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Abstract 

Additive white noise may significantly increase the response of bistable systems to a 
periodic driving signal. We consider two classes of double-well potentials, symmetric 
and asymmetric, modulated periodically in time with period 1/e, where e is a moder- 
ately (not exponentially) small parameter. We show that the response of the system 
changes drastically when the noise intensity a crosses a threshold value. Below the 
threshold, paths are concentrated near one potential well, and have an exponentially 
small probability to jump to the other well. Above the threshold, transitions between 
the wells occur with probability exponentially close to 1/2 in the symmetric case, and 
exponentially close to 1 in the asymmetric case. The transition zones are localised in 
time near the points of minimal barrier height. We give a mathematically rigorous 
description of the behaviour of individual paths, which allows us, in particular, to 
determine the power-law dependence of the critical noise intensity on e and on the 
minimal barrier height, as well as the asymptotics of the transition and non-transition 
probabilities. 

Date. December 29, 2000. 

2000 Mathematics Subject Classification. 37H99, 60H10 (primary), 34F05, 34E15 (secondary). 
Keywords and phrases. Stochastic resonance, noise-induced synchronization, double-well poten- 
tial, additive noise, random dynamical systems, non-autonomous stochastic differential equations, 
singular perturbations, pathwise description, concentration of measure. 

1 Introduction 

Since its introduction as a model for the periodic appearance of the ice ages [BPSV], 



stochastic resonance has been observed in a large number of physical and biological sys- 
tems, including lasers, electronic circuits and the sensory system of crayfish (for reviews 
of applications, see for instance | MW|1 ). 



The mechanism of stochastic resonance can be illustrated in a simple model. Consider 
the overdamped motion of a particle in a double-well potential. The two potential wells 
describe two macroscopically different states of the unperturbed system, for instance cold 
and warm climate. The particle is subject to two different kinds of perturbation: a 
deterministic periodic driving force (such as the periodic variation of insulation caused 
by the changing eccentricity of the earth's orbit), and an additive noise (modeling the 
random influence of the weather). Each of these two perturbations, taken by itself, does 
not produce any interesting dynamics (from the point of view of resonance). Indeed, the 
periodic driving is assumed to have too small an amplitude to allow for any transitions 
between the potential wells in the absence of noise. On the other hand, without periodic 
forcing, additive noise will cause the particle to jump from one potential well to the other 
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at random times. The expected time between transitions is given asymptotically, in the 
small noise limit, by Kramers' time, which is proportional to the exponential of the barrier 
height H over the noise intensity squared, namely e H ^ a . When both perturbations are 
combined, however, and their amplitudes suitably tuned, the particle will flip back and 
forth between the wells in a close to periodic way. Thus the internal noise can significantly 
enhance the weak external periodic forcing, by producing large amplitude oscillations of 
the system, hence the name of resonance. 

The choice of the term "resonance" has been questioned, as "it would be more appro- 
priate to refer to noise-induced signal-to-noise ratio enhancement" [Fox]. In the regime 
of a periodic driving whose amplitude is not a small parameter, one also speaks of noise- 
induced synchronization [SNAS|. Appreciable, though still sub-threshold amplitudes of 
the periodic driving have the advantage to enable transitions for small noise intensities, 
without requiring astronomically long driving periods. 

While the heuristic mechanism of stochastic resonance is rather well understood, a 
complete mathematical description is still lacking, though important progress has been 
made in several limiting cases. Depending on the regime one is interested in, several 
approaches have been used to describe the phenomenon quantitatively. The simplest ones 
use a discretization of either time or space. When the potential is considered as piecewise 
constant in time, the generator of the autonomous case can be used to give a complete 
solution [BPSV], showing that resonance occurs when driving period and Kramers' time 
are equal. Alternatively, space can be discretized in order to obtain a two-state model, 
which is described by a Markovian jump process |ET|. The two-state model has also been 
realised experimentally by an electronic circuit, called the Schmitt trigger | FH , McNWfl , 

In physical experiments, one has often access to indirect characteristics of the dynamics, 
such as the power spectrum, which displays a peak at the driving frequency. The strength 
of the resonance is quantified through the signal-to- noise ratio (SNR), which is proportional 
to the area under the peak (this definition obviously leaves some liberty of choice). The 
SNR has been estimated, in the limit of small driving amplitude, by using spectral theory 
of the Fokker-Planck equation [Fox, JH], or a "rate" equation for the probability density 
[McNW]. The signal-to-noise ratio is found to behave like e~ H ^ a /a 4 , which reaches a 
maximum for a 2 = H/2. 

The probability density of the process, however, only gives part of the picture, and 
a more detailed understanding of the behaviour of individual paths is desirable. Some 
interesting progress in this direction is found in [[ftj . The approach applies to a very general 
class of dynamical systems, in the limit of vanishing noise intensity. When the period of 
the forcing scales like Kramers' time, solutions of the stochastic differential equation are 
shown to converge to periodic functions in the following sense: The L p -distance between 
the paths and the periodic limiting function converges to zero in probability as the noise 
intensity goes to zero. Due to its generality, however, this approach does not give any 
information on the rate of convergence of typical paths to the periodic function, nor 
does it estimate the probability of atypical paths. Also, since the period of the forcing 
must scale like Kramers' time, the assumed small noise intensity goes hand in hand with 
exponentially long waiting times between interwell transitions. 

In the present work, we provide a more detailed description of the individual paths' 
behaviour, for small but finite noise intensities and driving frequencies. We consider two 
classes of one-dimensional double-well potentials, symmetric and asymmetric ones. The 
height of the potential barrier is assumed to become small periodically, which allows us to 
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Figure 1. A typical solution of the SDE (1.1) in the case of the symmetric potential (1.2) 



Heavy curves indicate the position of the potential wells, which approach each other at 
integer times. The straight line is the location of the saddle. Parameter values e = 0.01, 
a = 0.08 and ao = 0.02 belong to the regime where the transition probability between 
wells is close to 1/2. We show that transitions are concentrated in regions of order *Jo 
around the instants of minimal barrier height. 



consider situations where the period need not be exponentially large in 1/cr 2 for transitions 
between the wells to be likely. 

In the case of an asymmetric potential, we are interested, in particular, in determining 
the optimal noise intensity as a function of the driving frequency and the minimal barrier 
height, guaranteeing a close-to-periodic oscillation between both wells. We will estimate 
both the deviation (in space and time) of typical paths from the limiting periodic function, 
and the asymptotics of the probability of exceptional paths. The case of a symmetric 
potential shows an additional feature. For the right choice of the noise intensity, transitions 
become likely once per period, at which time the "new" well is chosen at random. We will 
again estimate the deviation from a suitable reference process and the asymptotics of the 
probability of exceptional paths. 

The systems are described by stochastic differential equations (SDEs) of the form 

d 

dx s = -— V{x s ,s)ds + adW s , (1.1) 

where W s is a Brownian motion. The potential V(x,s) is 1/e-periodic in s, and admits 
two minima for every value of s. The frequency e, the minimal barrier height between 
the wells and the noise intensity a are considered as (moderately) small parameters, the 
relation between which will determine the transition probability. 

The first class of potentials we consider is symmetric in x. A typical representative of 
this class is the potential 

1 1 

V(x, s) = —-a(es)x 2 + -x 4 1 , with a(es) = an + 1 — cos(2-7res). (1-2) 

Here an ^ is a parameter controlling the minimal barrier height. We introduce the slow 
time t = es for convenience. The potential has two wells, located at ±-\/ a(t), separated by 
a barrier of height |a(i) 2 . The distance between the wells and the barrier height become 
small simultaneously, at integer values of t. 
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Our results for symmetric potentials can be summarized as follows: 

In the deterministic case a = 0, we describe the dependence of solutions on t, ao and 
e (Theorem |2.1| ). Solutions starting at x > are attracted by the potential well at 
\J a(t), which they track with a small lag. If ao ^ e 2 / 3 , this lag is at most of the order 
e/ao; if ao ^ e 2 / 3 , it is at most of the order e 1 / 3 , but solutions never approach the 
saddle closer than a distance of order e 1 / 3 (even if ao = 0). 

When noise is present, but a is small compared to the maximum of ao and e 2 / 3 , 
the paths are likely to track the solution of the corresponding deterministic differ- 
ential equation at a distance of order <r/max{|t|, ^/oq, e 1 / 3 } (Theorem |2.2| ). The 
probability to reach the saddle during one time period is exponentially small in 
<7 2 /(max{a ,£ 2//3 }) 2 . 

If a is larger than both ao and e 2 / 3 , transitions between potential wells become likely, 
but are concentrated on the time interval [—y/a, y/a] (repeated periodically). During 
this time interval, the paths may jump back and forth frequently between both poten- 
tial wells, and they have a typical spreading of the order a/ maxjy^, e 1//3 }. After time 
y/a, the paths are likely to choose one of the wells and stay there till the next period 
(Theorem ^4|). The probability to choose either potential well is exponentially close 
to 1/2, with an exponent of order <r 3 / 2 /e, which is independent of ao (Theorem 2.3). 



• This picture remains true when a is larger than both yfaQ and e 1 / 3 , but note that 
the spreading of paths during the transition may become very large. Thus increasing 
noise levels will gradually blur the periodic signal. 

These results show a rather sharp transition to take place at a = maxjao, e 2 ^ 3 }, from a 
regime where the paths are unlikely to switch from one potential well to the other one, to 
a regime where they do switch with a probability exponentially close to 1/2 (Fig. |l]). 
The second class of potentials we consider is asymmetric, a typical representative being 

V(x,s) = -^x 2 + ^x 4 - X(es)x. (1.3) 

This is a double-well potential if and only if |A| < A c = 2/(3\/3). We thus choose X(es) = 
X(t) of the form 

X(t) = -(A c - a ) cos(2vrt). (1.4) 

Near t = 0, the right-hand potential well approaches the saddle at a distance of order 
ao, and the barrier height is of order a . A similar encounter between the left-hand 



potential well and the saddle occurs at t = 1/2. 

Our results for asymmetric potentials can be summarized as follows: 

• In the deterministic case, solutions track the potential wells at a distance at most of 
order min{e/ao, \/s}. If ao sS £, they never approach the saddle closer than a distance 
of order yfe (Theorem 2.5). 

• When a is small compared to the maximum of a 3 / 4 and e 3 / 4 , paths are likely to 

track the deterministic solutions at a distance of order a / max{y^[, aj , e x / 4 } (Theo- 

2//_r 3/4 



rem 2J: ) . The probability to overcome the barrier is exponentially small in a 2 /(max{a 3/4 ,e 3 / 4 }) 



For larger a, transitions become probable during the time interval [-a 2/3 ,cr 2/3 ]. Due 
to the asymmetry, the probability to jump from the less deep potential well to the 
deeper one is exponentially close to one, with an exponent of order <r 4//3 /e, while paths 
are unlikely to come back (Theorem |2.7| ) . 
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Figure 2 . A typical solution of the SDE (1.1) in the case of the asymmetric potential ( F3 ) 



The upper and lower heavy curves indicate the position of the potential wells, while the 
middle curve is the location of the saddle. Parameter values e = 0.005, a — 0.08 and 
ao = 0.005 belong to the regime where the transition probability between wells is close to 
1 . We show that transitions are concentrated in regions of order er 2 / 3 around the instants 
of minimal barrier height. 

• This picture remains true when a is larger than both a^ 4 and £ 1//4 , but the spreading 

of paths during the transition may become very large. 
Again, we find a rather sharp transition to take place, this time at a = max {a^ 4 , e 3 / 4 }. In 
contrast to the symmetric case, for large a the paths are likely to jump from one potential 
well to the other at every half-period (Fig. |2|) . 

In both the symmetric and the asymmetric case, we thus obtain a high switching 
probability between the potential wells even for small noise intensities, provided minimal 
barrier height and driving frequency are sufficiently small. They only need, however, to 
be smaller than a power of a: a® a and e <C cr 3 / 2 in the symmetric case, and ao <C a 4 / 3 , 
e <C ct 4 / 3 in the asymmetric case are sufficient conditions for switching dynamics. 

Our results require a precise understanding of dynamical effects, and the subtle in- 
terplay between the probability to reach the potential barrier, the time needed for such 
excursions, and the total number of excursions with a chance of success. In this respect, 
they provide a substantial progress compared to the "quasistatic" approach, which con- 
siders potentials that are piecewise constant in time. Note that some of our results may 
come as a surprise. In particular, neither the width (in time) of the transition zone nor 
the asymptotics of the transition probability depend on the minimal barrier height ao- In 
fact, the picture is independent of ao as soon as ao is smaller than e 2 / 3 (in the symmetric 
case) or e (in the asymmetric case), even for ao = 0. This is due to the fact that when 
ao is small, the time during which the potential barrier is low is too short to contribute 
significantly to the transition probability. 

The remainder of this paper is organized as follows. The results are formulated in detail 



in Section Subsection 2.2 being devoted to symmetric potentials, and Subsection 2.3 
to asymmetric potentials. Section ||| contains the proofs for the symmetric case, while 
Section ||] contains the proofs for the asymmetric case. 
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2 Results 

2.1 Preliminaries 

We consider non-autonomous SDEs of the form (|l.l| ). Introducing the slow time t = es 
allows to study the system on a time interval of order one. When substituting t for es, 
Brownian motion is rescaled and we obtain an SDE 

1 CT 

dx t = -f(x t ,t)dt + —=dW t , x to =x , (2.1) 
e ye 

where / is the force, derived from the potential V , and {W~t}t^t is a standard Wiener 
process on some probability space (Q, T , P). Initial conditions xq are always assumed to 
be square-integrable with respect to P and independent of {Wt}t^t - Without further 
mentioning we always assume that / satisfies the usual (local) Lipschitz and bounded- 
growth conditions which guarantee existence and pathwise uniqueness of a strong solution 
{xt]t of ( |2.1| ). Under these conditions, there exists a continuous version of {xt]t- Therefore 
we may assume that the paths uj \—* xt(uj) are continuous for P-almost all u G J7. 

We introduce the notation P* '^ for the law of the process {xt}t^t , starting in xo 
at time to, and use E* '^ to denote expectations with respect to F to ' x °. Note that the 
stochastic process {xt}t^t Q is an inhomogeneous Markov process. We are interested in first 
exit times of xt from space-time sets. Let icKx [to, t\] be Borel-measurable. Assuming 
that A contains (xo,io) ; we define the first exit time of (xt,t) from A by 

r A = M{t G [t ,ti]: (x t ,t)$A}, (2.2) 

and agree to set t^(uj) = oo for those uj G £1 which satisfy (xt(uj), t) G A for all t G [to, t\]. 
For convenience, we shall call T4 the first exit time of xt from A. Typically, we will 
consider sets of the form A = {(x,t) G R x [io>ii] : 9i(t) < x < #2^)} with continuous 
functions g\ < <?2- Note that in this case, T4 is a stopping timeQ with respect to the 
canonical filtration of (f2, J 7 , P) generated by {xt}t^t Q - 

Before turning to the precise statements of our results, let us introduce some notations. 
We shall use 

• \y\ for y ^ to denote the smallest integer which is greater than or equal to y, and 

• yV z and y A z to denote the maximum or minimum, respectively, of two real numbers 
y and z. 

• If ip{t, e) and ip{t, e) are defined for small e and for t in a given interval /, we write 
ij)(t,e) x (p(t, e) if there exist strictly positive constants c± such that c~Lp{t,e) ^ 
vp(t, e) ?J c+(p(t, e) for all t G / and all sufficiently small e. The constants c± are 
understood to be independent of t and e (and hence also independent of quantities 
like a and ao which we consider as functions of e). 

• By g(u) = 0{u) we indicate that there exist 5 > and K > such that g(u) ^ Ku 
for all u G [0, 5], where 5 and K of course do not depend on e or on the other small 
parameters ao and a. Similarly, g(u) = o(l) is to be understood as ]im u ->o g(u) = 0. 

Finally, let us point out that most estimates hold for small enough e only, and often only 
for P-almost all w£Jl. We will stress these facts only where confusion might arise. 

1 For a general Borel-measurable set A, the first exit time ta is still a stopping time with respect to the 
canonical filtration, completed by the null sets. 
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2.2 Symmetric case 



We consider in this subsection the SDE (2.1) in the case of / being periodic in t, odd in 
x, and admitting two stable equilibrium branches, with a "barrier" between the branches 
becoming small once during every time period. A typical example of such a function is 



f{x,t) = a(t)x 



with a(t) = ao + 1 — cos 2irt. (2.3) 
£ , which we assume to satisfy 



We will consider a more general class of functions / : R 2 - 
the following hypotheses: 

• Smoothness: f G C 4 (.M,M), where Ai = [—d,d] x R and d > is a constant; 

• Periodicity: f(x,t + 1) = f(x,t) for all (x, t) £ Ai; 

• Symmetry: f(x,t) = —f(—x,t) for all (x,t) £ A4; 

• Equilibrium branches: There exists a continuous function x* : R — ► (0, d] with the 
property that f(x,t) = in Ai if and only if x = or x = ±x*(i); 

• Stability: The origin is unstable and the equilibrium branches ±x* are stable, that is, 
for alii e R , 



a(t) :=d x f(0,t) > 
o*(t) :=d x f(x*(t),t) <0. 



(2.4) 



Behaviour near t = 0: We want the three equilibrium branches to come close at integer 
times. Given the symmetry of /, the natural assumption is that we have an "avoided 
pitchfork bifurcation", that is, 



cW(0,0) < 

a{t) = a + ait 2 + 0(t 3 ) 



(2.5) 



where a± > and d xxx f{0, 0) are fixed (of order one), while ao = ao(e) = o e (l) is a 
positive small parameter. Is is easy to show that x*(t) behaves like \/a(t) for small t, 
and admits a quadratic minimum at a time t* = O(ao). Moreover, a*(t) x — a(t) near 
t = 0. 

We can choose a constant T £ (0,1/2) such that the derivatives of a(t) and x*(t) 
vanish only once in the interval [— T, T]. We finally require that x*(t), a(t) and a*(t) 
are bounded away from zero outside this interval. We can summarize these properties 
as 



x*(t) 

a(t) 
«*(*) 

We start by considering the deterministic equation 

j„det 

e^ r = f(xr,t). 



' Vdo 


for \t\ ^ y^Oo 




< \t\ 


for ^ \t\ 


^ r 


1 

V 


for T < t ^ 1 


— T 


ao 


for t ^ y^Oo 




< t 2 


for ^ \t\ 


^ r 


1 

V 


for T < t < 1 


— T 


-o(t) 


for all t. 





(2.6) 

(2.7) 
(2.8) 

(2.9) 
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Without loss of generality, we may assume that xf et starts at time — 1+T in some x^ +T > 
0. Tihonov's theorem [|Grj , [Ti| , applied on the interval [—1 + T, — T], implies that xf et 
converges exponentially fast to a neighbourhood of order e of x*(t). We may thus assume 
that x^, = x*(—T) + 0(e). In fact, since x* is decreasing at time — T, we may even 
assume that x_£, — x*(—T) X e. 

The motion of xf et in the interval [— T, T] is described in the following theorem. 

Theorem 2.1 (Deterministic case). The solution xf et and the curve x*(t) cross once 
and only once during the time interval [-T, T] . This crossing occurs at a time t satisfying 
i—t*x (e/ao) A e 1 / 3 . There exists a constant cq > such that 



„dct 



X*(t) 



J for-T^t ^-co^VeVS) 



(2.10) 



and thus xf et X |i| in i/tese iime intervals. For \t\ ^ co(^/ao V e 1 ^ 3 ), 



„det 



a i/ a ^ e 2/3 
x / 3 tfa ^e 2 / 3 . 



(2.11) 



Finally, the linearization of f at xf et satisfies 

a(t) := d x f(xf ct ,t) X -{t 2 V a V e 2 / 3 ). (2.12) 



We give the proof in Subsection |3.1| . The relation ( |2.11 ) may be surprising, since it 



means that no matter how small we make ao, xf et never approaches the saddle at x = 
closer than a distance of order e 1//3 . This fact can be intuitively understood as follows. 
Even if ao = and near t = 0, we have 



j det r det 

e—— > -const (x det ) 3 =► ^f et > const *° (2.13) 

(,/ ^l + (x^) 2 (t-t )/e 



Since x^ ct X e 1 / 3 for to ~ — e 1 / 3 , xt cannot approach the origin significantly during any 
time interval of order e 1 / 3 . After such a time, however, the repulsion of the saddle will 
make itself felt again, preventing the solution from further approaching the origin. In 
other words, the time interval during which a(t) is smaller than e 2 / 3 is too short to allow 
the deterministic solution to come close to the saddle. 

We return now to the SDE (|2.1| ) with a > 0. Assume that we start at some deter- 
ministic x_i + t > 0. Theorem 2.3 in [ |BG| ] shows that the paths are likely to track the 
deterministic solution xf et with the same initial condition at a distance of order cr 1-5 for 
any 5 > (with probability ^ 1 — (1/e 2 ) exp{— const/a 26 }), as long as the equilibrium 
branches are well separated, that is, at least for — 1 + T ^ t ^ —T. A transition between 
the potential wells is thus unlikely if a = 0(|loge|~ 1//25 ), and interesting phenomena can 
only be expected between the times — T and T. Upon completion of one time period, i. e., 
at time T, the Markov property allows to repeat the above argument. Hence there is no 



limitation in considering the SDE (2.1) on the time interval [— T, T], with a fixed initial 
condition x-t satisfying x^t — x*(—T) x e. We will denote by xf et and xt, respectively, 
the solutions of (|2.9D and (|2.1|) with the same initial condition x_t- 
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FIGURE 3. Solutions of the SDE (2.1) with symmetric drift term (2.3), shown for two 




different noise intensities, but for the same realization of Brownian motion. Heavy curves 
represent the equilibrium branches ±x*(i), and the straight l ine represents the saddle. 
Smooth light curves are solutions of the deterministic equation ( |2.9| ) tracking the potential 
wells with a small lag, while rugged curves are paths of the SDE. Parameter values are 
e = 0.01, a = 0.02, a = 0.02 (left) and a = 0.08 (right). 

Let us start by describing the dynamics in a neighbourhood of xf et . The main idea is 
that for a sufficiently small, the typical spreading of paths around xf et should be related 
to the variance v(t) of the solution of Q2.1| ), linearized around xf et . This variance is given 
by 

a 2 r* - r* 

v (t) = — e 2a(M)/£ ds, where a(t, s) = / a{u) du. (2.14) 

£ J-T Js 

The variance is equal to zero at time — T, but behaves asymptotically like a 2 /\2a(t)\. In 
fact, if we define the function 

C(t) := ^~ T)/£ +- f ^ S)l£ ds, (2.15) 

2\a{-T)\ eJ_ T 

then v(t) differs from (7 2 ^(t) by a term that becomes negligible as soon as \a(t, — T)\ is 
larger than a constant times e|loge|. C(t) has the advantage to be bounded away from 
zero for all i, which avoids certain technical problems in the proofs. We shall show that 

r V ao V £ z ' 6 

We introduce the set 



B(h) = {(x,t): \t\ ^ T,\x-xf et \ < h^/Oj)}, (2.17) 
and denote by tjsqa the first exit time of xt from B(h). 

Theorem 2.2 (Motion near the stable equilibrium branches). There exists a con- 
stant Hq, depending only on f, such that 

• if -T < t < -(y/aS V e 1/3 ) and h < h t 2 , then 



-T,x. 



{r m <t}^C(t,s)expL~ i-0(e)-o(±\ }; (2.18) 
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*/ -{y/ao V e 1 / 3 ) < f ^ T and h < h (a V £ 2/3 ) ; then 



-T,x. 



In both cases, 



f 1 h 2 

\j~B(h) < t} < C(i,e)exp<^ ~2^2 



1 - 0(e) - O 



Of) 



Ve 2 / 3 



(2.19) 



C(t,e) 



-T)| +2. 



(2.20) 



We give the proof in Subsect ion |3.2j . This result has several consequences. Observe first 
that the exponential factors in fl2.18|) and (|2.19|) are very small as soon as h is significantly 
larger than a. The prefactor C(t,e) (which, unlike the exponent, we do not believe to 
be optimal) leads to sub exponential corrections, which are negligible as soon as h/cr > 
C(|loge|). It mainly accounts for the fact that the probability for a path to leave B(h) 
increases slowly with time. The theorem shows that the typical spreading of paths around 



xf et is of order 



a 



a 



\t\ V Joq~ 



V eV3 ' 



(2.21) 



If a <C ao V £ 2//3 , we may choose h S> cr for all times, and thus the probability of leaving a 
neighbourhood of xf et , let alone approach the other stable branch, is exponentially small 
(in er 2 /(ao V e 2 / 3 ) 2 ). On the other hand, if a is not so small, ( p. 18 ) can still be applied to 
show that a transition is unlikely to occur before a time of order —\/a. Figure ^ illustrates 
this phenomenon by showing typical paths for two different noise intensities. 

Let us now assume that a is sufficiently large to allow for a transition, and examine 
the transition regime in more detail. We will proceed in two steps. First we will estimate 
the probability of not reaching the saddle at x = during a time interval [to>*l]- The 
symmetry of / implies that for any t ^ t\ and xq > 0, 



¥ to ' Xo {x t < 0} 



1 



' x °{3s 6 (to, t) : x s = 0} 



2 

= * _ l^ X0 {x s >0VsG [t ,t]} 
> I_I P ^o{ Xs >0VsG [t ,ti]}. 



(2.22) 



In the second step, we will show, independently, that paths are likely to leave a neigh- 
bourhood of x = after time \fa. Thus if the probability of not reaching x = is small, 
the probability of making a transition from the positive well to the negative one will be 
close to 1/2 (it can never exceed 1/2 because of the symmetry). This does not exclude, 
of course, that paths frequently switch back and forth between the two potential wells 
during the time interval [—y/a, y/a]. But it shows that ( 2.22 ) can indeed be interpreted 
as a lower bound on the transition probability. 
Let 5 > be a constant such that 



xd xx f(x,t) < 



for Ixl ^ 5 and \t\ ^ T. 



(2.23) 



Our hypotheses on / imply that such a 5 of order one always exists. In some special cases, 
for instance if f(x,t) = a{t)x — x 3 , S may be chosen arbitrarily large. 
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Theorem 2.3 (Transition regime). Let c\ > be a constant and assume c\a ^ ao V 
e 2 / 3 . Choose times —T ^ to ^ *i ^ T with t\ G [— cxy/a, c\yfa\, and let h > 2a be 
such that Xg et + h\J £(s) < 5 for all s G [to,ii]- Then, if c\ is sufficiently small and 
x e(0,xt^ + ^h^/ctto}], 

F t0 ' X0 {x s > Vs G [to,h]} 

+ A e -*fc"V +2 exp{-K - * ^i,-ci^) | (2>24) 
2 V e / I log(/i/<7) e J 

where R is a positive constant, and a(t, s) = j a(u) du. 

The proof is given in Subsection 3.3. The first term in ( |2.24 ) is an upper bound on the 
probability that xt escapes "upward". Indeed, our hypotheses on / do not exclude that 
other stable equilibria exist for sufficiently large x, which might trap escaping trajectories. 
The second term bounds the probability of x s remaining between and x^ et + h^J £(s) for 
—c%y/a ^ s ^ t\. This estimate lies at the core of our argument, and can be understood 
as follows. Assume x s starts near x^ et . It will perform a certain number of excursions 
to attempt reaching the saddle at x = 0. Each excursion requires a typical time of order 
As, such that a(s + As, s) >c e (that is, a(s)As x e), in the sense that the probability of 
reaching before time s+As is small. After an unsuccessful excursion, x s may exceed if et , 
but will return typically after another time of order As. Thus the total number of trials 
during the time interval [— c\yfo, t\] is of order a(ti, —ciy/a)/e. Under the hypotheses of 
the theorem, the probability of not reaching the saddle during one excursion is of order 
one, and thus the total number of trials determines the exponent in ( 2.24 ). 

Before discussing the choice of the parameters giving an optimal bound in Theorem |2.3| . 
let us first state the announced second step, namely the claim that the paths are likely to 
escape from the saddle after t = c\\fo. For k G (0, 1), let us introduce the set 

V{k) = |(x,t) G [-5, 5} x [ci V^, T] : > na{t) } . (2.25) 

The upper boundary of V(k) is a function x(t) = \/l — k (1 — 0(t))x*(t). Let Tx>t K ) denote 
the first exit time of xt from D(k). 

Theorem 2.4 (Escape from the saddle). Let < k < 1 and assume c\a ^ oq V e 2 / 3 . 
Then there exist constants C2 ^ c\ and Cq > such that 

/ t \2 e -Ka{t,t 2 )/2e 

P t2 '* 2 {^M > *} < c Arl , WTW > (2 ' 26) 

W \\Jo-J y?l _ e -K.a(t,t 2 )/e 

for all (x2,t2) G T>(k) with t<i ^ c<i\fo ■ 



The proof is adapted from the proof of the similar Theorem 2.9 in BG |. Compared 
to that result, we have sacrificed a factor 2 in the exponent, in order to get a weaker 
condition on a. We discuss the changes in the proof in Subsection 3.4 . 

For the moment, let us consider t2 = C2\fo. We want to choose a t such that a(t, ^2) ^ 
e|logcj|. Since a(t, t2) is larger than a constant times i 2 (i — £2), it suffices to choose a t of 
order \/ct(1 + e|log o"|/cr 3//2 ) for ( |2.26| ) to become small. Hence, after waiting for a time of 
that order, we find 

^ t2,X2 {rv( K ) > t} < const |loga|W 2 , (2.27) 
which shows that most trajectories will have left T>(k) by time t, see Fig. BL 
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FIGURE 4. A typical path xt of the SDE (2.1) in the symmetric case, shown in a neigh- 
bourhood of the origin, where the potential barrier reaches its minimal height. We show 
here a situation where the noise intensity is large enough to allow a transition. The po- 
tential wells at ±x*(t) behave like ±(y/ao V |i|). The deterministic solution xf et starting 
near the right-hand potential well tracks x*(t) at a distance at most of order (s/ao) V e 1 / 3 , 
and never approaches the saddle at x — closer than ^/ao Ve 1 / 3 . The path xt is likely 
to stay in the set B(h) up to time —y/h when h 3> a. Between times — \fo and y/a, the 
path is likely to reach the origin. It may continue to jump back and forth between the 
potential wells up to time y/a, but is likely to leave a neighbourhood T> of the saddle for 
times slightly larger than yfa. For each realization w such that xt(uS) reaches the saddle 
at a time r, there is a realization u> such that xt(u>') = —xt(to) is the mirror image of Xt 
for t ^ r, which explains why the probability to choose one well or the other after the 
transition region is close to 1/2. 



It remains to show that the paths are likely to approach either x*(t) or — x*(t) af- 
ter leaving T>(k). Let us first consider the solution x t det of the deterministic differential 
equation ( |2.9D with initial time ti ^ ci\fo and initial condition \x2 — x*(t)\ cti for 
some small constant c > 0. Here we need to choose c small in order to arrange for 
a(t) = d x f(xf et ,t) X — t 2 which allows us to proceed as in our investigation of the mo- 
tion for t ^ — (y/a V e 1 / 3 ), cf. Subsections 3.1 and 3.2. Under these assumptions, xf et 
approaches a neighbourhood of ±x*(t) exponentially fast and then tracks the equilibrium 
branch at distance e/t 2 . As before, one can show that the path xt of the solution of the 
SDE ( |2.lD with the same initial condition is likely to remain in a strip around xf et of 
width scaling with a j \J"a{t) X cr/t. So if a path xt leaves T>(n) at time Tx>r K ), then this 
path is likely to approach x*(t), if x Tt}{k) is positive, and —x*(t), otherwise. Note that 
l x *( r x>(re)) — x TT3t K ) I nas to be smaller than CTp( K ), which restricts the possible values for k. 
Therefore, we choose k small enough to guarantee that x*(t) — x(t) ^ ct for all t ^ c\\fo. 
Finally note that paths which are not in T>(k) at time ci\fa but are not further away from 
=tx* than ct at some time t will also approach the corresponding equilibrium branch. 
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Let us now discuss the choice of the parameters in ( 2,24 ) giving an optimal bound. 



For a ^ (ao V e 2 / 3 )/c 2 , Theorem 2.2 shows that up to time t slightly less than —yfo, the 



paths are concentrated around xf et . Therefore, here we should choose an initial time to 
slightly before —\fo. In order for the second term in ( p. 24 ) to be small, we want to choose 



t\ — (—c\y/a) as large as possible. Note, however, that h^J £(s) has to be smaller than 
5 — Xg 6t for all s. Since the order of ((s) is increasing for s ^ — (y/a V e 1 / 3 ), it turns out 
to be more advantageous to choose t\ negative and of order —y/a, say t\ = —\c\^fo. In 
this case, \a(lu\,—c\\fG*)\ is larger than a constant times cr 3//2 (independently of ao), and 
this does not improve significantly for larger admissible t\. At the same time, this choice 
allows us to take h x 5y/a. We find 



{x s > Vs G [t ,ciy/a]j < P*°' x °|x s > Vs G [to, -\c x ^a}\ (2.28) 



^, NV^e-^+exp/ C0USt ^ 



a J e \ \og{5 2 /a) e 

Consider first the generic case <5xl. The second term in ( |2.28| ) becomes small as soon as 
c/(|logcr|) 2 / 3 » e 2//3 holds in addition to the general condition c\a ^ ao V e 2 ^ . The first 
term is small as long as a = 0(1/ log(cr/,s 2 / 3 )). We thus obtain the following regimes: 

• for a ^ ao V e 2 / 3 , the transition probability is exponentially small in a 2 /(ao V e 2//3 ) 2 ; 

• for a ^ ao/c 2 with (e|loge|) 2 / 3 <C a l/|loge|, the probability of a transition between 
the wells is exponentially close to 1/2, with an exponent given essentially (up to 
logarithmic corrections) by 

a 3 / 2 1 

— A-; (2.29) 

£ <J 

• for a ^ l/|loge|, the paths become so poorly localized that it is no longer meaningful 
to speak of a transition probability. 

( 2.29 ) shows that the transition probability becomes optimal for a x e 2 / 5 . For larger 



values of the noise intensity, the possibility of paths escaping "upward" becomes sufficiently 
important to decrease the transition probability. However, if the function / is such that 5 



can be chosen arbitrarily large, the second term in ( 2.29 ) can be removed without changing 
the first one (up to logarithmic corrections) by taking 5 2 = cr/e, for instance. In that case, 
transitions between the wells become the more likely the larger the ratio cr/e 2 / 3 is. 

One should note that a typical path will reach maximal values of the order a^J C(0) x 
o'Av^o V e 1 / 3 ). Thus, due to the flatness of the potential near t = 0, if a is larger 
than y/ao V e 1 / 3 , the spreading of the paths during the transition interval is larger than 
the maximal distance between the wells away from the transition. In general we cannot 
exclude that paths escape to other attractors, if the potential has more than two wells. 

It may be surprising that the order of the transition probability is independent of ao 
as soon as a > ao- Intuitively, one would rather expect this probability to depend on 
the ratio a 2 ,/ a 2 , because of Kramers' law. The fact that this is not the case illustrates 
the necessity of a good understanding of dynamical effects (as opposed to a quasistatic 
picture). Although the potential barrier is smallest between the times — ^/oq and ^/oq, 
the paths have more opportunities to reach the saddle during larger time intervals. The 
optimal time interval turns out to have a length of the order ^/a, which corresponds to 
the regime where diffusive behaviour prevails over the influence of the drift. 
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2.3 Asymmetric case 

We consider in this subsection the SDE Q2.1| ) in the case of / being periodic in t and ad- 
mitting two stable equilibrium branches, but without the symmetry assumption. Instead, 
we want each of the potential wells to approach the saddle once in every time period, but 
at different times for the left-hand and the right-hand potential well. A typical example 
of such a function is 

f(x,t) = x - x 3 + X(t) with X(t) = -(A c - a )cos2vrt. (2.30) 

Here A c = 2/(3\/3) is defined by the fact that / has two stable equilibria if and only if 
|A| < A c . Observe that d x f vanishes at x = ±x c = ±l/vo, and 



f(x c + y, t) = A c - (A c - ao) cos 2vrt - V3y 2 - y 3 

V 



a + 27r 2 (Ac-a )t 2 + C'(t 4 )-v / 3y 2 "' 3 



Here the function f(x c ,t) plays the role that a(t) played in the symmetric case, and near 
t = the right-hand potential well and the saddle behave like x c ± 3 _1 / 4 y / /(x c , t), while 
the left-hand potential well is isolated. Near t = 1/2, a similar close encounter takes place 
between the saddle and the left-hand potential well. 

We will consider a more general class of functions / : R 2 —* R , which we assume to 
satisfy the following hypotheses: 

• Smoothness: f G C 3 (.M,R), where A4 = [—d, d] x R and d > is a constant; 

• Periodicity: f(x,t + 1) = f(x,t) for all (x,t) G Ai; 

• Equilibrium branches: There exist continuous functions x*_ < Xq < x+ from R to 
[— d, d] with the property that f(x,t) = in M. if and only if x = x±(i) or x = x^t); 
the zeroes of / should be isolated in the following sense: for every 8 > 0, there 
should exist a constant p > such that, if \x — ^ 5 and \x — Xq(*)I ^ ^> then 

/(•'••/! ->/'.Q 

• Stability: The equilibrium branches are stable and the equilibrium branch Xq is 
unstable, that is, for all t € R , 

a* ± (t):=d x f(x* ± (t),t)<0 

a* o (t):=d x f(x*(t),t)>0. [ ' ' 

• Behaviour near t = 0: We want x\ and Xq to come close at integer times. Here the 
natural assumption is that we have an "avoided saddle-node bifurcation", that is, 
there exists an x c £ (— 5, 5) such that 

d xx f{x c ,0)<0 
8 x f(x c ,t) = 0(t 2 ) (2.33) 
f{x c ,t) = a + ai t 2 + 0(t 3 ), 

where a\ > and d xx f(x c ,0) are fixed (of order one), while ao = ao(e) = c e (l) 
is a positive small parameter. These assumptions imply that x+(t) reaches a local 
minimum at a time t* + = O(oq), and Xq(£) reaches a local maximum at a possibly 



2 Since / depends on a small parameter ao, we want to avoid that f(x,t) approaches zero elsewhere but 
near the three equilibrium branches, even when ao becomes small. 
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different time tg = 0{clq). We can assume that for a sufficiently small constant T > 0, 
the three equilibrium branches and the linearization of / around them satisfy 

x* (t) - x - { ^ f ° r ^ a* (t) x /~V^o f ° r I*' ^ ^ 

C \ \t\ for ^ ^ |i| ^ T + \-\t\ for ^/oJi < |t| < T 

x*(t) - x x J~^° fOT 1*1 ^ ^ *m x J ^ for 1*1 ^ ^ 

C \-\t\ for^^|i|^T ° \ \t\ for ^Hq" < |i| < T 

x* (t) - z c x -1 for |t| < T a* (t) x -1 for \t\ < T. (2.34) 

• Behaviour near t = t c : We want x* and Xq to come close at some time t c £ (T, 1 — T). 
This is achieved by assuming that similar relations as ( |2,33j ), but with opposite signs, 
hold at a point (x' c ,t c ). 

• Behaviour between the close encounters: To exclude the possibility of other almost- 
bifurcations, we require that x* (i) — Xg(t) and Xg(i) — x* (i), as well as the derivatives 
( |2.32j ), are bounded away from zero for T < t < t c — T and t c + T<t<l — T. 

Note that a sufficient assumption for the requirements on the behaviour near (x' c ,t c ) to 
hold is that f(x, t + |) = —f(—x, t) for all (x, t). 

We start by considering the deterministic equation 

j det 

e^ r = f{xr,t). (2.35) 

As in the symmetric case, it is sufficient to consider the dynamics in the time interval 
[— T, T], with an initial condition satisfying x^ — x*^(—T) x e. The situation in the time 
interval [t c — T,t c + T] can be described in exactly the same way. 

Theorem 2.5 (Deterministic case). The solution xf et and the curve x+(t) cross once 
and only once during the time interval [— T, T] . This crossing occurs at a time t satisfying 
i — t+x (e/y/ao) A \fe. There exists a constant cq > such that 



„dct 



<(*) 



j for-T^t^-c {y/a^V^) 
-||r forc {^V ^t^T, 



(2.36) 



and thus xf et — x c x \t\ in these time intervals. For \t\ ^ co(y / ao V \fe), 

4 et -x c x(^ 2/a °^ £ (2.37) 
[y/e ifa ^e. 

The linearization of f at xf ct satisfies 



a(t) :=d x f(x? c \t) x V Vao V v^)- (2-38) 

Moreover, (|2.35| ) admits a particular solution x^ et tracking the unstable equilibrium branch 
Xg(t). It satisfies analogous relations, namely, x t det and Xg(t) cross once at a time i 
satisfying i- t% x -(t - i^), and (g3p , (^37|) and (|2.38|) /ioW for x det and xj$(i), but 
with opposite signs. 
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The proof is similar to the proof of Theorem 2.1 , and we comment on a few minor 
differences in Subsection 4.1. Note that ( |2 .37 ) implies that xf et never approaches the 
saddle at Xg(i) closer than a distance of order yfe. 

We return now to the SDE fl2.1| ) with a > 0. We will denote by xf et and xt, respectively, 
the solutions of ( [2.35 ) and (|2.l|) with the same initial condition x-t satisfying x-t — 
x*,(—T) X e. We introduce again the function 



2|o( 



1 e m(t,-T)/e + 1 /"* e 2a(M)/ e dSj where g(tj s) = /" fi( 

— )l £ J-T Js 



u)du, (2.39) 



which behaves, in this case, like 

C(t) - n i= j= for 1*1 < T - (2.40) 

We define once more the set 

B(h) = {{x,t): \t\ < T,\x-xf et \ < /i^CW}, (2.41) 
and denote by TjgrM the first exit time of from B{h). 

Theorem 2.6 (Motion near the stable equilibrium branches). There exists a con- 
stant ho, depending only on f, such that 
• */ — T < t < — (y^oo V y/e) and /i < h \t\ 3 / 2 , then 



~ T ' x - T {r B (h) <t}< C(t,e)expi-~ 



if -(V^O V Ve) ^t s^T and h < h Q {al ,4: V e 3 / 4 ), tfien 



-T,x_ 



1 /l 2 

'{TH(ft) < *} < C(t,e)exp<j -- — 



1 - 0(e) - O 



3/4 



Ve 3 / 4 



In both cases, 



C(t,e) = -\a(t, -T)|+2. 



(2.42) 



(2.43) 



(2.44) 



This result is proved in exactly the same way as Theorem It has similar con- 
sequences, only with different values of the exponents. The typical spreading of paths 
around xf et is of order 



a 



(2.45) 



If a <C a^ 4 Ve 3 / 4 , the probability of leaving a neighbourhood of xf et , or making a transition 
to the other stable equilibrium branch, is exponentially small (in a 2 /(a^ 2 Ve 3 / 2 )). On 
the other hand, if a is not so small, ( ^42| ) can still be applied to show that a transition is 
unlikely to occur before a time of order — <r 2 / 3 . 

Let us now assume that a is sufficiently large for a transition to take place, i.e. that 
a ^ a^ 4 V e 3 / 4 . We want to give an upper bound on the probability not to make a 
transition. Let us introduce levels 5q < 8% < x c < 62 such that 



d xx f{x,t) ^ 



for 5o ^ x ^ 01 and \t\ ^ T 
for Si ^ x ^ 82 and \t\ ^ T. 



(2.46) 
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FIGURE 5. A typical path xt of the SDE (2.1) in the asymmetric case, shown near t = 0, 
where the right-hand well at x* + approaches the saddle at Xq. We show again a situation 
where the noise intensity is large enough to allow a transition. The deterministic solution 
xf et starting near the right-hand potential well tracks x+(t) at a distance at most of order 
(e/y/ao) A y/e, and never approaches the saddle at Xq closer than ^/ao V \fe. The path x t 
is likely to stay in the set B(h) up to time —h 2 / 3 when ft> a. Between times —ti 2 / 3 and 
h 2 / 3 , the path is likely to reach the saddle. Although it may fall back into the right-hand 
potential well, it is likely to finally overcome the potential barrier and reach a level Si of 
order 1 below the saddle, after which it quickly reaches a lower level 5q. The distance 
between 5\ and <5o can be much larger than in this picture. Finally, the path will track the 
deterministic solution starting in x = 5o, which approaches a neighbourhood of order e of 
the left-hand potential well at x*_. 



Here 5q and 5\ are always of order 1 (in fact, we must have 8q > x* (£) for all t £ 
[— T, T]), and we think of 8o as being in the basin of attraction of x*_. Our hypotheses 
imply that a 82 of order one satisfying (2.46) always exists, but 82 may be chosen arbitrarily 
large in particular cases such as f(x, t) = x — x 3 + A(t). 

The non-transition probability can be estimated by distinguishing three cases: 

• Either xt, starting in xo > x c at some to < 0, never reaches 8\. The probability of this 
event can be shown to be small in a similar way as in Theorem |2.3| , the main difference 
being that due to the asymmetry, we can do better than estimating the probability 
not to reach the saddle. 

• If xt reaches 81 , one can estimate in a very simple way the probability not to reach 8q 
as well, using the fact that the drift term is bounded away from zero. 

• If xt reaches 8q, Theorem 2.3 in |BG| shows that xt is likely to reach a small neigh- 
bourhood of x*_ as well. 
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Theorem 2.7 (Transition regime). Letci andci he positive constants and assume that 



anc 



c{ a ^ a f Ve 3 / 4 . Choose times -T ^ to ^ i x < t < T with t x G [-cicr 2 / 3 , ci<7 2 / 3 ] 
t ^ ti + C2S. Let h > 2a be such that x^ et + /iv/ C(s) < #2 / or s £ [io>*i]- Then, for 
sufficiently small c\, sufficiently large C2 and a// xo £ (£i, xf et + |^-\/ C(*o)]> 

P^{x. > 5 V, G [to,t]} £ * + 1) e-i*-/- 3 

3 r 1 S(ti,-cio- 2 / 3 )-i 

+ 2 eXP l" K log(/ l /^) V|loga| e J 

+ e- K/,j2 , (2.47) 

where k is a positive constant, and a(t,s) = J"* d x f(x^ et , u) du. 

The proof is given in Subsection The three terms on the right-hand side of ( p. 47 ) 
bound, respectively, the probability that x t escapes through the upper boundary x^ ct + 
h*J C(s), the probability that xt reaches neither the upper boundary nor 6%, and the 
probability that xt does not reach So when starting on 5\ (Fig. |5|). The crucial term is the 
second one. 

Let us now discuss the optimal choice of parameters. If we choose t\ = — tjCict 2 / 3 , we 
can take h >c ^a 1 ^ 3 , where 62 = d~2 — x c , and we get the estimate 



P t0 < X0 {x s ><5 Vs G [t 0> *]} 



^ £ e -o(^) + / ^1) + e -*,o> . (2 . 48) 

e I log(5 2 /<7 4 / 3 ) V|loga| e J V ; 

As in the symmetric case, when 52 X 1, we obtain the following regimes: 

• for <7 ^ Oq^ 4 \/£ 3 / 4 , the transition probability is exponentially small in a 2 / (og^Ve 3 / 4 ) 2 ; 

• for ag^ 4 V e 3 / 4 <C cr <C (l/|log e|) 3 / 4 , the transition probability is exponentially close 
to 1, with an exponent given essentially (up to logarithmic corrections) by 

(T 4 / 3 1 

• for <t ^ l/|loge|, the paths become so poorly localized that it is no longer meaningful 
to speak of a transition probability 

The transition probability becomes optimal for a X e 3//§ . Note, once again, that the 
exponent is independent of ao- 

If the function / is such that 62 can be chosen arbitrarily large, the second term in 
( 2.49 ) can be removed without changing the first one (up to logarithmic corrections) by 
taking <5>| = cr 8//3 /e for instance. If a > a^ 4 V £ 1//4 , the paths may become extremely 
delocalised in the transition zone, and could escape to other attractors. 
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3 Symmetric case 



We consider in this section the nonlinear SDE 

dx t 



-f(xt,t)dt + -^=dW t , 



(3.1) 



where / satisfies the hypotheses given at the beginning of Subsection 2.2. By rescaling x, 
we can arrange for d xxx f(0, 0) = —6, so that Taylor's formula allows us to write 



where go , g\ £ C 3 satisfy 



f(x,t) = x[a(t) +g (x,t)) 
d x f{x,t) = a(t) + gi(x,t), 



9o(x,t) = [-l + r (x,t) 
gi(x,t) = [-3 + n(x,t) 



(3.2) 



(3.3) 



with continuously differentiable functions ro,r% satisfying ro(0, 0) = r\(0, 0) = 0. 

The implicit function theorem shows the existence, for small t, of an equilibrium curve 



x*(t) = [l + 0{^Mt))]^/a{t). 



(3.4) 



For small t, the curve x*(i) behaves like ^/oo [l + 0((t/ \/ao) 2 ) + 0(- s /ao) l , and it admits a 
quadratic minimum at some time t* = 0{clq). Thus we can choose a constant T £ (0, 1/2) 
such that 



x*(t) 



^Eq for \t\ ^ ^Jao 
\t\ for y/ao ^ \t\ ^ T 
1 for T < t< 1 - T. 



(3.5) 



3.1 Deterministic case 

In this subsection we consider the deterministic equation 

s^ = f(x u t). (3.6) 
at 

As already mentioned, Tihonov's theorem allows us to restrict the analysis to the time 
interval [—T,T], and to assume that X-t — x*(—T) >c e. 

Remark 3.1. During the time interval [— T, T], the process xt crosses the equilibrium 
branch x*(t) once and only once, the time t of the crossing satisfying t ^ t*. This fact 
is due to the property that xt is strictly decreasing when lying above x*(t), and strictly 
increasing when lying below. Let t% < t2 < • • • be the times of the successive crossings of xt 
and x*(t) in [-T, T]. Then xt is decreasing between — T and t\ (since X-t — x*(—T) > 0), 
increasing for t\ < t < t2, and so on. Thus x*(t) must be increasing for t slightly larger than 
t\, and decreasing for t slightly larger than t^. Since, by assumption, x*(t) is decreasing 
on [— T, t*) and increasing on (i*, T], this implies that t\ ^ t* and ti > T. Therefore, there 
is at most one crossing. We shall see below that x% and x*(t) actually cross and we will 
also determine the order of that time t. 
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We consider now the difference yt = xt — x*(t). It satisfies the equation 



dy j., , , da;* 

£-^ = a*(t)y + &*(y,t)-e— 



where Taylor's formula, ( |3.2| ) and ( |3.3| ) yield the relations 
a*(t) 



&*(y,*) 



dx* , . 



-2a(t)[l + 0(Va(^)] x 
-(3x*(t) + y)y 2 [l + 0(x*(t) + y)] 



ao for |t| ^ -^ao 
t 2 for y^o < |*| < T 



-1 for -T < t < -^/oq 

(t- t* 



1 



for |i| ^ 
for yfoQ ^ t ^ T, 



(3.7) 

(3.8) 
(3.9) 

(3.10) 



with t* = t*(ao) = O(a ). 

We start by giving a technical result that we will need several times. 

Lemma 3.2. Let a(t) be a continuous function satisfying a(t) x — ((3 V t 2 ) for \t\ ^ T, 
where (3 = (3(e) ^ 0. Let xo x 1> ar ^ define a(t,s) = a(u) du. Then 



/J V e 2 / 3 
1 

f2 



/or |i| < y^Ve 1 / 3 



/or^Ve 1 ' 3 < |t| < T. 



(3.11) 



Proof: To prove the lemma, we take advantage of the fact that the expression on the left- 
hand side of ( 3.1 1| ) is the solution of an ordinary differential equation. By the semi-group 
property, we may consider separately the following regimes: For f3 e 2 / 3 , we distinguish 
the cases t € [-T,-T/2], [-T/2,-^], [-^,-s/a ], [-e/a , e/a ], [e/a ,VP], [VP,T\ 
and for (3 < e 2 / 3 , we deal separately with t £ [-T, -T/2], [-T/2, -e l, \ [-e 1 / 3 ,-^], 
[-\/^,V?], [v^iS 1 ^ 3 ]) [s 1 ^ 3 ,T]. On each of these time intervals the claimed behaviour 
follows easily by elementary calculus, see also the similar result [BG, Lemma 4.2]. □ 



Proposition 3.3. There exists a constant cq > such that the solution of (|3.7| ) with 
initial condition y_T x e satisfies 

y t x^ for-T^t^ -c (^V£ 1/3 ). (3.12) 

Proof: Let Co ^ 1 and c\ > T 2 y-r/£ be constants to be chosen later, and denote by r 
the first exit time of yt from the strip < yt < cie/t 2 . Set to = —cq{\/<H) V e 1 / 3 ). Then, 
for -T < t < r A t , we get from (O) and flO) that 



|&*(y,i)| 3x*(t) + y y , / e \ e . / cA c\ 



yt 2 |t| |tp V A |*IV l<l 3 V Co/cg 
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for some constants M,M' > 0. The relations ( [3.8] ) and ( 3.10| ) yield the existence of 
constants c± > such that a*(t) ^ — c_t 2 and — -£x*(t) ^ c + for t G [— T, to]. From (3.7) 
and fl3~13l) we obtain 



d 2/ / .2 

£ dt ^ " C - t * 



l-^fl + - C ^ 



+ ec+. (3.14) 



°0/ °0- 

For any given ci, we can choose Co large enough for the term in brackets to be larger than 
1/2. Then, by Lemma 3^, there exists a constant C2 = C2(c + , c_) > such that 

Vt < y_ T e- c -(* 3 + T3 )/ 6e +c + y%- c -( t3 - s3 )/ fe d S < c 2 ^ (3.15) 

for all t £ [— T, r A to]- Therefore, if c\ > C2, then r ^ to follows. 

The lower bound can be obtained in exactly the same way. □ 

For the remainder of this subsection, let to = — cq{^/clq Ve 1 ' 3 ) with cq chosen according 
to the preceding proposition. Note that this proposition implies that xt x x*(t) x |t| for 
— T ^ t ^ to and, in particular, that yt >c (e/ao) A e 1 ^ 3 . 

We now consider the dynamics for |t| ^ \to\, starting with the case of ao not too small, 
i. e., the case of yt Q x e/«o- 

Proposition 3.4. There exists a constant 70 > 0, depending only on f andyt , such that, 
when ao ^ 7o£ 2//3 , then 

yt = C 1 (t)(t* -t) + C 2 (t) with Ci(t)x-^, C 2 (t)x-^ (3-16) 
/or a// |t| ^ |to|. 

Proof: Again, we will only show how to obtain an upper bound, since the corresponding 
lower bound can be established in exactly the same way. 

First we fix a constant c\ > aoyt Q /e + 2(t* — to)/y/ao + 4e/(c_a ). We denote by r 
the first exit time of yt from the strip \yt\ < c\eja§. For to ^ t ^ r A |to|, we have 



\b*(y,t)\ 



^M(3x*(t) + \y\)\y\ < M'a ( 1 + c™) c™ (3.17) 



|y| 

with constants M, M' > 0. Choosing 70 large enough, we get 



dv C— s 
£ -£4 t -—a y-c--={t-t*), (3.18) 
at 2 Ja$ 



which implies 



ft 

yt < y* e~ c - ao( *-* o)/2£ / e - c - ao( *- s)/2E (s - 1*) ds 

x/OQ ,/t 



^L( t * _ t ) + ^( e ) e -c_ao(*-*o)/2 £ + (3 19) 
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by partial integration. Here 

77(e) 



Vto 



2^ 
a 



to 



+ 



2s 



C-CL. 



3/2 




(3.20) 



satisfies 77(e) = 0(e/ao). 

We want to estimate the contribution of the middle term on the right-hand side 
of fl3,19p . Assume first that 77(e) > and consider t ^ i*. By convexity, 



-c~ao{t—to)/2£ 



t*-t 



Now, 



Since X e~ x -^OasI 



2e 

3/2 



+ 77(e) 



to 
1 



+ e" 



-c-ao(t*—to)/2e 



t* 



to 



3/2 ' 



■ 00, we also have 



+ - 



4e 2 



5/2 




a, 



5/2' 




(3.21) 



(3.22) 



(3.23) 



provided 70 is large enough. This shows the existence of constants C\ ^ 2 and C2 > 
such that 

- e - e 2 

y* ^Ci — (i*-i) + C 2 — ^ for t^t* At. (3.24) 



For t > t* 



3/2 



-c_a (i-to)/2e 



5/2 



^ e 



-c_a (i*-to)/2e 



(3.25) 



is immediate, and (3.23) shows that (3.24) also holds for t* ^ t ^ r. Note that in the 
case 7/(e) ^ 0, ( 3.24 ) holds trivially. Since yt < c\e/ao is a direct consequence of ( |3.19D 
and our choice of c%, r ^ |to| follows, and, therefore, the upper bound (|3T2l ) holds for all 
|t| ^ |i |. □ 



Note that the result ( 3.16 ) implies that yt changes sign at a time t* + O{e/ao), which 
shows that xt actually crosses x*(t) at a time t satisfying t — t*x e/clq. For large enough 
70, the proposition also shows that xt x y/ao for \t\ ^ |io| and that y\ to \ x — e/ao- 

We consider now the case ao < 7o£ 2 ^ 3 with 70 from Proposition 3.4. Without loss of 
generality, we may assume that 70 ^ 1. 

Proposition 3.5. Assume that ao < 7oe 2//3 . Then, for any fixed t\ x e 1 / 3 , 

,1/3 



x t x £*'" /or t < i ^ *i, 
and xt crosses x*(t) at a time t satisfying t x e 1 / 3 . 

Proof: In order to show ( 3.26j ), we rescale space and time in the following way: 



(3.26) 



1/3 1/6 ; 



x = e ' a/ z 



t = e 



V3„-l/3 ( 



Let sq = e 



-i/3„,y 3 



a x to- Then z so 
dz 



1, and z satisfies the differential equation 



ds 



(s,e)z+ [-l + ro( 



eVy/« Zje VV 1/8 



(3.27) 



(3.28) 
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where 

a(s, e) = \ /3 = a + , 2 + 0{e^s\ a = — ^ < . (3.29) 

( 3.28 ) is a perturbation of order e 1 / 3 of the Bernoulli equation 

dz 

— = a(s)z — z 3 , a(s) = a$ + s 2 . (3.30) 
as 

Using Gronwall's inequality, one easily shows that on an s-time scale of order 1, the 
solution of (|3.28| ) differs by 0(e 1 ^ 3 ) from the solution of ( |3.30D , which is 



7 p a(s,so) 



1/2 ' 

l + 2zl I e 2 ^"' s °W 
so 



where a(s, sq) = / a{u)du. (3.31) 

J sn 



This function is bounded away from zero, and remains of order one for s of order one, 
which shows that xt x e 1 ^ on [in,ii]. 

Since xt x s 1 ^ 3 and x*(t) X y/aQ V \t\, xt and x*(i) necessarily cross at some time 
tx^. □ 

Note that the above proposition also implies bounds on yt, namely, yt = 0{e 1 / 3 ) for 
to ^ t K Hi an d there exist constants c+ > c_ > such that 

' eV3 for t ^ i< c_e 1 /3 ) 
^ x < for t = i, (3.32) 
^-e 1 / 3 for c+e 1 / 3 < t ^ ii. 

Gathering the results for ao ^ 7o£ 2 ^ 3 an d ao < 7o £2//3 i we see that there exists a time 
ti x (y^oVe 1 / 3 ) such that y tl x — e/if. By enlarging en if necessary, we may assume that 

*1 =Q)(V^V£ 1 / 3 ). 



Proposition 3.6. On £/ie interval [t\,T\, 



y t x (3.33) 



£ 

Proof: The proof is similar to the one of Proposition |3.3| . □ 

Note that the previous result implies Xt x x*(t) x t, provided Co is large enough. 
So far, we have proved that for t G [—T,T], x t tracks x*{t) at a distance of order 

^ A — Ae 1 / 3 , (3.34) 
t z ao 

and that the two curves cross at a time i satisfying t — t* x (e/ao) A e 1 / 3 . Let us now 
examine the behaviour of the linearization 

a(t)=d x f(x t ,t), (3.35) 

which will determine the behaviour of orbits starting close to the particular solution xt- 
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Proposition 3.7. For all t £ [—T,T] and all ao = o £ (l), 

a(t) x -(t 2 Vo Ve 2/3 ). (3.36) 
Proof: By Taylor's formula we get 

a{t) =d x f(x*(t)+y t ,t) 

= a*(t) - [6 + 0(x*(t) + y t )] (ar*(t) + ^t)^- (3.37) 

Consider first the case ao ^ 7o£ 2 ^ 3 for 70 large enough. Equati on (|3.8| ) implies that 
a*(t) — c_(ao V i 2 ) for a constant c_ > 0. On the other hand, ( |3.34D shows that the 
second term on the right-hand side of ( |3.37| ) is bounded in absolute value by c+e/ao for a 
constant c + > 0. Thus if 70 > (c + /c_) 2 / 3 we obtain that a(t) X a*(t) X — (ao V t 2 ). 

We consider next the case ao < 7o£ 2 ^ 3 - For \t\ ^ co£ 1//3 , the above argument can 
be repeated. The non-trivial case occurs for \t\ < cq£ 1//3 . By rescaling variables as in 



Proposition 3.5, we obtain that 

a(t) = e 2/3 a\ /3 [a{s, e) - 3z 2 + 0(e 1/3 )] . (3.38) 

We have to show that a(t) x — e 2 / 3 which is equivalent to a(s, e) — 3z 2 x —1 for s of order 
one. The lower bound is trivial as d(s,e) ^ and z s x 1. In order to show the upper 
bound, first note that for t ^ 0, we have x% > x*(t) which implies a(t) < a*(t). Therefore, 
it is sufficient to consider s ^ 0. Taking into account the expression (|3.3l| ), we find that 
showing the upper bound amounts to showing that 

const + ^) e - 2fi(s ' So) +2d(s) f e~ 2&(s ' u) du < 3. (3.39) 

Since |s | is proportional to cq, choosing a priori a large enough cq also makes |so| large. 
Thus it is in fact sufficient to verify that 

2a(s) f e" 2 ^'") du < 3 (3.40) 
J —00 

for all s ^ 0. Optimizing the left-hand side with respect to ao ^ and s shows that we 
may assume do = and that ( p. 40 ) holds. □ 



3.2 The random motion near the stable equilibrium branches 

We now consider the SDE 

dx t = -f(x t ,t)dt+-^=dW t , X- T = xo, (3.41) 

e ye 

on the time interval [— T, T], where we assume xq — x*(—T) x e. In order to compare the 
solution xt with the solution x^ et of the corresponding deterministic equation ( |3.6| ), we 
introduce the difference y t = x t — xf et , which satisfies the SDE 




(3.42) 
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where a(t) is the linearization ( |3.35| ) of / along xf et , and Taylor's formula yields the 
relations 



b(y, t) = -[l + 0(xf ct + \y\) + 0(t)] (3xf ct + y)y 2 , 
\b(y,t)\^M(xt ct + \y\)y 2 



(3.43) 



whenever \t\ ^ T and xf ct + \y\ ^ d, where M is a positive constant. Let us first consider 
the linearization of (|3.42|), namely 



dy t ° = -a{t)y Q t dt + ^ dW t , y°_ T = 0. 
e ye 



(3.44) 



The random variable y® is Gaussian with expectation zero and variance 



v(t) 



— I e 2a(M)/£ ds, where a(t, s) = f 

£ J-T Js 



a(u) du. 



Lemma 3.2 and Proposition fO] imply that 

1 



C(t) ■■- 



1 



e 2a(t,-T)/e + 



1 



e 2ffi(t,a)/e x — — . 

e £ 2 Vo V e 2 / 3 



(3.45) 



(3.46) 



2|o(-T)| 

Thus u(t) is of order cr 2 /(t 2 V ao V £ 2//3 ), except for t very close to — T. We now show that 
y° is likely to remain in a strip of width proportional to y 7 C(i). 

Proposition 3.8. For -T < t < T, 

,01 , , lfc 2 



-T,0 



sup 



I//,' 



^ /i} < C(i, e) exp{---2 (1 - 0(e)) }, (3.47) 



where 



C(t, e ) = M^ + 2. 



(3.48) 



Proof: Let — T = no < ui < • • • < uk = with some K > 0, be a partition of [— T, i]. 
In [BG, Lemma 3.2], we show that the probability ( 3.47j ) is bounded above by 



K ( lh 2 

2 Pk, where P k = exp\ inf 

k—i 



coo 

C(«fe) 



3 2a(u fc ,s)/e 



}. (3.49) 



Now we choose the partition by requiring that 



a(u k ,u k - 



-2e z 



for 1< < if 



\<*(t,-T)\ 
2e 2 



Since a(s) < 0, we have ('(s) = [2a(s)£(s) + l]/e ^ 1/e, and thus 



inf «£U 1 



inf 



tifc - S 



s 



(3.50) 



(3.51) 
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If k is such that \uk\ y/<M)V e 1 ^ 3 , then by ( |3. 50 ) and ( 3.36 ), there is a constant c_ such 
that 



2e z ^ c / s 2 ds ^ — u 2 k (u k - u k -i), 



and hence by ( 3.46j ) (choosing the same c_ for brevity of notation) 

u k - u k -i I2e 2 u\ 

5^- 



For all other fc, we have 

2e 2 > c_(a \/e 2/3 )(u k - u fc _i) 
In both cases, we find 

P fc ^exp{-^(1-C»(e))} ; 
which leads to the result, using the definition of K. 



u k - n fc _i 2£ 
C(itjfc) ^ c 2 



(3.52) 

(3.53) 

(3.54) 

(3.55) 
□ 



Let us now compare solutions of the linear equation (3.44) and the nonlinear equation 



( 3.42 ). We introduce the events 



Q t (h) = {u: \y s \ ^ hyXis) Vs G [-T,t]} 
n° t (h) = {u: \y° s \ < h^/Os) Vs G [-T,t]}. 



(3.56) 
(3.57) 



Notation 3.9. For two events Q\ and Q.2, we write J7i C 0,2 if ^-almost all u> G fii 
belong to f^- 

Proposition 3.10. There exists a constant g, depending only on f , such that 
• if-T^t^ -(y/aS V e 1/3 ) and h < t 2 /g, then 

h 



1 + £> 



i 2 



; 



if-iy/ao Ve 1/3 ) s^t^T andh< (a V e 2/3 )/g, then 

h 



SVi(h) c fi< 



i + e- 



(3.58) 



(3.59) 



a Ve 2 / 3 - 

Proof: The proof is based on the fact that the variable z s = y s — y° s satisfies the relation 



1 

e J-t 



(3.60) 



Consider first the case — T ^ t ^ — (y/ao V e 1 / 3 ). Let p > be a constant to be chosen 
later, and set 5 = gh/t 2 < 1. We define the first exit time 



r = inf{s G [-T,t]: \z s \ ^ 5hy/((sj} G [-T,t] U {oo}. 



(3.61) 
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Pick any u G A := Qf(h) n {uj : t{uj) < 00} and s G [— T, t(w)]. Then we have 

I^MI^/tTcR, |y„(o;)|< (1 + 5)/i\/cR < 2/i^cR (3.62) 

for all u G [— T, s]. From ( |3.5| ) and ( |3.46| ), we obtain the existence of a constant c+ > 
such that 



xt ct < c+|n|, < ( 3 .63) 



for these it. Hence, by ( 3.43j ) we get the estimate 

\/c+\,,?c+ ,,h 2 c 2 + f 2h 
4/i 2 4 <4Af-i-±(l + -— 
|u| / u z |s| V y / C+S 

and thus, by ( 3.60| ) and Lemma [3^ , 

\ Zs \<4M ^4fl + _^-)l f e 5 ( s ' u >/ e d U <4M^(l + — 

where we use again the same c+ for brevity of notation. Using ( |3.46| ) once again, we arrive 
at the bound 

<4Af-^=^(l + —->). (3.66) 



|6(l/«,«)| < M(c + |-u| +2/i^-f )4/i 2 -t ^4M^-±( 1 + ^=3) (3.64) 



Now we choose 



which implies 



2 c 3 

V 8M — — , (3.67) 



q h f h\ 5 



< ^{ 1 + ^)^2 {1 + S)<5 (3 - 68) 



for all s G [— T, r(u;)] ; by the definition of 5. Hence |« T ( W )| < Shy/ ((t(uj)) for almost all uj G 
A. Since we have |-z r (u)l = Shy/ £(t(u)) whenever t(u) < 00, we conclude that ¥(A) = 0, 
and thus t(uj) = 00 for almost all uj G Q®(h), which implies that |y s (w)| ^ (1 + S)hy/ ((s) 
for — T ^ s ^ t and these uj. This completes the proof of ( 3,5§| ), 

The proof of ( |3.59|) is almost the same. In the case — (y/ao V e 1//3 ) ^ t ^ T, we take 
(5 = gh/(ao V e 2 / 3 ). The estimate ( |3.63| ) has to be replaced by 



xt ct ^ c+(\u\ V^Ve 1 / 3 ), \y u (uj)\ < 2fe v _2 (3.69) 

|it| V y'ao V e 1 ' 6 

and thus we get, instead of ( |3.65| ), the bound 

h 2 c% ( 2h \ . . 

\z s \ ^4ikf — - — - ±- l + _____ _ (3.70) 



/^VeV3)( s 2 V a Ve 2 /3) ^ ^(a Ve 2 / 3 ). 
The remainder of the proof is similar. □ 



Now, the preceding two propositions immediately imply Theorem as Proposi- 
tion shows the desired behaviour for the approximation by a Gaussian process and 



Proposition 3.1C allows to extend this result to the original process. 
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3.3 The transition regime 

We consider now the regime of a sufficiently large to allow for transitions from one stable 
equilibrium branch to the other. Here xf et is the solution of the deterministic equation 
Q3.6D with the same initial condition as in the previous sections, which tracks x*(t) 
at distance at most 0(e 1//3 ). xt denotes a general solution of the SDE (|3.1|). Our aim is 
to establish an upper bound for the probability of not reaching the axis x = 0, which, 
by using symmetry, will allow us to estimate the transition probability. Let 5 > be the 
constant defined in ( 2.23 ), i.e. by 

xd xx f(x,t) < for \x\ 5 and \t\ < T. (3.71) 

The basic ingredient of our estimate is the following comparison lemma which allows 
us to linearize the stochastic differential equations under consideration and, therefore, to 
investigate Gaussian approximations to our processes. The lemma gives conditions under 
which relations between initial conditions carry over to the sample paths. 

Lemma 3.11. Fix some initial time to G [ — T,T]. We consider the following processes 
on [t ,T]: 

• the solution xf et of the deterministic differential equation (3J3) with initial condition 
< ct e[0,<5]; 

• the solution xt of the SDE ( 3.4l| ) with initial condition xt G [a;^,^]; 

• the difference yt = %t~ xf et , which satisfies yt = Xt — xf et 0; 

• the solution y®of the linear SDE 

dy° t = -a(t)y® dt + dW t , where a(t) = d x f(xf et ,t) (3.72) 
e ye 

with initial condition y® G [yt ,<5 — xf^]. 
If ^ y° ^ 5 — x^ ct for all s G [to,t], then y s ^ y® for those s. Similarly, if ^ y s ^ 
5 — Xg et for all s G [to,t], then y Q s ^ y s for those s. The result remains true when t is 
replaced by a stopping time. 



Proof: The hypothesis fl3.71| ) implies that for all y £ [0,5 — r 



,detl 

y *z \y,u — a> 

det i „. „\ ^ jv^det 



f{xT + y,s)^ f(xf\ s) + a(s)y. (3.73) 



Let r = inf{s G [t ,t]: y s g" [0,5 - xf ct ]} G [t ,t] U {oo}. For t s ^ r, the variable 



z s = ys — satisfies 



*h + - f [f(4 et + Vu,u)- f(xt et ,u) - a(u)y° u ] du 

e J to 



'to 
1 f s 

< Zt + - / a(u)z u du, z to ^O. (3.74) 
e Jto 

Applying Gronwall's inequality, we obtain 

z s < zto e &{s ' to)/e ^ Vs G [to, t At], (3.75) 

where a(s, to) = J* t ^ a{u) du. This proves the result for to ^ s ^ r At. Now if y s is negative, 
the result is trivially satisfied, and if y s becomes positive again, the above argument can 
be repeated. Note that y T ^ 5 — x^ et is immediate. This proves the first assertion, and 
the second assertion can be proved directly, without use of r. □ 
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We will now proceed as follows. Let ((t) be the function defined in ( |3.46 ), and let h 
be such that xf et + h-\J Q(s) < 5 for all s G [to,t]. Given xq G (0,6), we can write 

„det 



F to ' xo {x s > Vs G [to,*]} < ^ to ' xo \ sup Xs J!l > h\ 



+ P'°' XO {0 < x s < x** + h^/C(s) Vs G [to,t]}. (3.76) 

We will estimate these two terms separately. The first event is similar to the event we have 
examined in the previous subsection, but we need here an estimate valid for all times, even 
when a is not very small, whereas the previous result is only useful for a ^ t 2 V ciq V e 2 / 3 . 
We will show the following. 

Proposition 3.12. Assume s$ x < xf Q et + \h^/C(tQ). Then 



dot 



sup 



5 / 1«(*, * )| 



+ 1 e 



(3.77) 



where k is a positive constant and a(t,to) = J* a(s) ds. 



Proof: 

1. We define a partition to = uq < ui < • • ■ < uk = t of the interval [to, £] by requiring 

|«(*, *o)| 



for 1 < k < K 



Note that similar arguments as in the proof of Proposition |3.q yield 

Uk+l ~ Uk 



C(Uk) 



0(e) 



for all k. 



(3.78) 



(3.79) 



Now let pk = \h\J C( u k) an d y s = x s — Xg et as usual. Define 

9Uk,Vk) ■■■■■ V s 



Qk = sup 

Vk^Pk 



for < k < K - 1, and 



{ 



sup 



> h 



} 



+ - 



VC(s) 

u k ,y k ( gup V^^^y >p \ 



Qk-i = sup 

Vk-i^Pk-i 



tUK-ltVK-l 



sup 



> h 



}■ 



(3.80) 



(3.81) 



Then 



dct 



sup 



} 



o{ sup -^_>fcj + p*o.Wo{ sup < /», y Ul > Pl | 



+ E t °'%>{l { y „ 1<Pl} P Ul '^{ sup 



>"}} 



K-l 
fc=0 



(3.82) 
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(k) 

2. In order to estimate Qk, we introduce the stochastic process (y s ) s &[u k u k+1 ] defined 
by " ~' ' ' 

y( k ) = p k e «( s > u k)/e + ° f S e a(s,u)/e dw (k) ^ (3^3) 
V £ Ju k 

where {W^) u ^[u k ,u k+1 } ls the Brownian motion = W u — W Uk . Note that y^> is 
the solution of the SDE ( |3.72| ) with initial condition y^ k = Pk at time Uk- We define 
the stopping times 

r° = inf{s e [tifc.Ufc+i]: yf ] = 0} 

. . (3.84) 

r+ = inf{s G [ufc,t* fc+ i] : j/W = h^({s) } 

describing the time when y^ either reaches the t-axis or the upper boundary h^J ((s). 

Now, Lemma p. 1 1| implies that if y Uk ^ pk, then y s ^ y s for Uk ^ s ^ r A r + . This 
shows that 



Q fc < p«*>A*{ T ° < Uk+l ) + P «*,P*{ T + < Ufc+1 } +P « fc ,P fc {y(fc) +i > pfc+1 } (3. 8 5) 



»«fe,Pfc| T ?( , , , \ _|_ p«fe,Pfc/T-+ ^ ,,, , , \ _i_ F"fc,Pfe/ g ,(fc) 

for ^ k < K — 1, and 

Q K ^ ^ p^-i>te-i| r o < UA -} +p«a-i,pk-i| t + < UK j ( 3 86 ) 

Each of these terms depends only on y( k \ and can be easily estimated. Let 

(g =— / e 2a(« fc+1 , U )/ edn ( 3> g 7 ) 



"A- 



denote the variance of yu k+1 - Then by symmetry (as in fl2.22| )), we have 

F Uk ' pk {r < u k+1 } = 2¥ Uk ' pk {yi k k l i < 0} 

u k+1 ,u k) /e( v (k) y 1 /* 

e- 2 / 2 dz 



2 r ^ k+1 ' Uk)/e ^y 112 , 



\p2/JX J-oo 

( \ p l e 2a(«*+i,«fc)/e 

^ exp ~2 — (*5 

* P l 8 6 ^H^' 



The second term on the right-hand side of ( |3.85|) or ( |3.86| ), respectively, can be esti- 
mated using the symmetry (in distribution) of ( p.83 ) under the map a i— > — a: 

^P«{r° <n fe+1 }. (3.89) 
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In order to estimate the third term on the right-hand side of ( 3.85 ), we will use the 
fact that for k < K — 1 



CK+i) = C(«fc) e 2a(uk +i> Uk)/£ +- / e Mu k +x,s)/e d 



1-e" 



inf 



2 «fc^s^M fc+1 |a(s)| 

Proposition [T?] and ( |3.46 ) thus yield the existence of a constant c_ > such that 



(3.90) 



C(«fc+i) ^ -2 , 1 - e " 



C(«fc) ' 2c 2 _ 

This allows us to estimate (for k < K — 1) 



(3.91) 



l(Pfe+i-/0fce 5 ^+i^)/ £ ) 2 



(fc) 



exp< 



h 2 ((u k )( /CK+i) e _i 



(fc) 



C(wfc) 



(3.92) 



3. The estimates ( ggD , fl3~89p and (ggD , inserted in (|3~85|) and fl3~86p , imply that 

5 r n 2 



o r ft i 
Qfc < 2 6XP l Kk a^J ' 



with 



1 CKO _2 



1 A , 1 



e 2 -l 



(3.93) 



(3.94) 



By ( 3.87 ), for each k, there exists a Ok G [e , 1] such that 



( fc ) c \ 

^V+i (itfc+i-ttfc) 



<7^ 



(3.95) 



Together with ( 3.79j ), this implies that x 1 for all fc, and thus the result follows 
from ( |3.82| ) with R = inf & . □ 

We now give an estimate of the second term in ( |3.76| ). The Markov property implies 
that we will obtain an upper bound by starting at time —c\\fa. 



Proposition 3.13. There exist constants c\ > and k > such that, if c\o ^ oq V e 2 ^ 3 
and h > 2a, then 



1 a(ti, -Cjy/a) \ 



^ 2 exp< —At 



log(ft/cr) 



7 



(3.96) 



ZioWs wtft — ci\fb~) = X^ Clv ^ a ( s ) ds, /or —c\\fo ^ ii ^ c\\fo and all initial condi- 
tions xq satisfying ^ xq ^ x l C c lv ^ + C( — civ^)- 
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Proof: 

1. Let q = g(h/a) ^ 1 and define a partition —c\y/a = uq < ■ ■ ■ < uk = ti of [—ciy/a,t\] 
by 



a(u k ,u k -i 



QE 



iorl^k<K 



a(ti, -ci^/a) 



Q£ 



(3.97) 



We would like to control the probability of not reaching the i-axis during the time 
interval [u k ,u k+ i]. Let 



Qk 



sup 



P tt *'**{0 < x s < x** + hy/C(8j Vs G [u k ,u k+1 ]}. (3.£ 

0<a !fc <xdrt + fty'C( Ufc ) 

Then the probability on the left-hand side of ( |3.96[) is 

— IB" —c\\fB ,a;o J i 

^ \ x {o<x s ^dct +h ^/^y vsel-civ^^K-i]} 

p «A--i,^ K _ 1 | 0<x ^ x det + /lv ^y Vs6[«jr_i,«jc]}} 

< Q A -_iP- Cl ^°{0 < x s < xf* + /iTcR Vs G [-d^^-i]} 



fsr-1 

< •• ■ < n 

fc=0 



(3.99) 



If we manage to estimate each Q k by a constant less than 1 (say, 1/2), then the 
probability will be exponentially small in K. In the sequel, we shall estimate Q k 
uniformly in k = 0, ... K — 2, and bound Qk-i by 1, since the last interval of the 
partition may be too small to get a good bound. So let k < K — 1 from now on. 
2. We consider first the case < x k ^ We define the process {xi) Uk ^.s^u k+1 as the 
solution of the linearized SDE 



dx^ fc) = a(s)4 fc) ds + 4= dWj fc) , 



(3.100) 



where (W< 



W l r 

(fe) 



is the Brownian motion Ws^ = W s — W Uk . Let v^ k+1 denote 



(k) 



the variance of Xu k+1 . Then 



-2a(ufc+i,u fc )/e (fc) 

1 



-2a 



(7 



inf 



1 



2 a(tt) 



[1 



-2o(n fc+ i,« fc )/el 



1 



-2 e 



<7 



(3.101) 



2 a (life) V a(tifc+i) ' 

We can now apply Lemma 3.11 in the particular = to show that if < 

x s ^ 5 for s G [ujfc, Mfc+i], then cc|^ ^ x s in the same interval. We thus obtain 

p«*,* fc { < x s < z?* + / lv / CM Vs G [u fc ,n fc+1 ]} < P«*.**{a;(*) > Vs G [u k ,u k+1 ]}. 

(3.102) 
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The probability on the right-hand side satisfies 

F u k , Xk { x (k) > Vs e ^ nfc+l] } = x _ 2P«*-**{4*2n < °) 

= 2P u *- a: *{x£) +i > 0} - 1, (3.103) 

yielding 



P™{0 < x s < * + hy/Os) Vs G 



-1/2 



2 



, Xke \ V ™k+l) y /? 

< -= / e~ z 12 dz 

2 x k 



/ p -2a(u k+1 ,u k )/e «,( fc ) 



2 / 1 \/fl(Ht)V«(%l) det 



By making ci small enough, we can guarantee that this bound is smaller than some 
imposed constant of order 1, say 1/2. This shows that the length of [u k ,u k+ i] has 
been chosen large enough that the probability of reaching the t-axis during this time 
interval is appreciable. 
3. We examine now the case x^ < x k < x~^ + hy/((s). We introduce a time u k G 
(u k ,u k+ i), defined by 

a(u k ,u k ) = -x0£- (3.105) 

Our strategy will be to show that xt is likely to cross x^ et before time u k , which will 
allow us to use the previous result. Proposition 3.7 implies the existence of a constant 
L > such that 

—Q£= / a(u)du^L / (—a(u))du = L\a(u k ,u k )\. (3.106) 

Ju k Ju k 

Let {y^)u k 4-s4-u k+1 be the solution of the linear SDE 



dyi k) = a(s)yW ds + ^ dW$ k \ yff =Vh = x k - x%, (3.107) 



(k)\ (k) 

where (W s ) s &[ Uk ,u k+1 ] is again the Brownian motion W$ = W s — W Uk . The variance 

r (k) . 
of yl k ' is 

fif = °- I'" Ptoxy* ds > 1 " e " g/L , , /, , (3.108) 
Uk e J Uk ' 2 \a(u k )\ V |o(« fc )| V ; 

Lemma 3.11 shows that if Xg Ct ^ x s ^ 5 on the interval [tijfc,Wfc], then x s — x^ et ^ 
on that interval. If we introduce the stopping time 

T k = infjs G [u k ,u k ]: x s = xf et } G [u k ,u k+1 ] U {oo}, (3.109) 
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then we have 



^ Xk {0<x s ^xf t + h^/C{s) V S £^,« H i]} 



<j p « fc ,x fc | x dct <Xg <: x det + h ^/Q(s) y s € [ Uk ,u k ]} 

+ E u ^{l {Tk<ilk} P Tk '< t {0 <x s ^ x^ + hy/^s) Vs e [r fc ,u fc+1 ]}}. (3.110) 
The second term on the right-hand side can be bounded, as in ( |3,104 h by 



1 



y/a(u k ) V a(u k +i) 



^/tt V l-es 



a 



sup x^ et . 

u k ^s^u k 



(3.111) 



Using ( |3.108| ), the first term on the right-hand side of ( 3. 110| ) can be estimated in the 
following way: 



iU k ,X k 



{xf l <x s ^ xf* + h^/CiT) Vs G [u h , %]} 

2 y ke a(u k ,u k )/e 



2tt fd^ft 



(C) 



h 



-\ 1 777\/|a(wfc)| V |a(ti fc )| -VC(^fc) e " 

■y/7r V 1 — e"" 1 cr 



-g/2L 



(3.112) 



Using Proposition and ( 3.46| ), it is easy to show that the expression (|a(«fc)| V 
\a(uk)\)Q{uk) is uniformly bounded by a constant independent of k and e. The sum of 
(3.111) and of the last term in (3.112) provides an upper bound for Q k - 
4. Using the fact that for |u| ^ c\\fo and c\a ^ aq V e 2 / 3 , one has a(u) = O(cla) and 



x 



det 



O(ciy/o), we arrive at the bound 



(3.113) 



where the constant C can be chosen independent of g because g ^ 1 by assumption. 
Thus if we choose cf ^ 1/4C and g = 2Llog(4Ch/a) V 1, we obtain that Qk ^ 1/2 for 
k = 0, . . . , K - 2. This yields 

< V < 2exp{-(log2) a(tl '- Cl ^ ) }, (3.114) 



fc=0 



and the result follows from our choice of g. 



□ 



Now the proof of Theorem 2.3 follows from ([3.76 ) and the two preceding propositions, 
where we use the Markov property to "restart" at time —c\^fa before applying Proposi- 



tion 3.13. 



3.4 Escape from the saddle 

In this subsection, we investigate the behaviour of the random motion xt given by the 
SDE ( [3.1| ) for t ^ t2 ciy/a, i.e., after the transition regime. We want to show that xt is 
likely to leave a suitably defined neighbourhood of the saddle within time C(e|log er| /£?,). 
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The proof of Theorem 24 is very similar to the proof of [BG, Theorem 2.9], and for the 
sake of brevity, we will refrain from giving all the details. Instead, we will discuss how to 
proceed and then focus on those parts which need to be modified. 

From now on, we will assume that t ^ t 2 ^ c\yfo and that a is large enough in order 
to allow for transitions, i. e., a ^ ao V e 2 / 3 . We want to estimate the first exit time Tx>( re ) 
of x t from the set 



V(k) = \ (x,t) G [- 



-6,5] x [ci^T]: > «a(t) 



(3.115) 



where k E (0, 1) is a constant. Note that the upper boundary x(t) of T>{n) satisfies 
x(t) = y/l — k(1 — 0(t))x*(t). Our first step towards estimating Td< k ) is to estimate the 
first exit time ts from a smaller strip S, defined by 



(x,t) G [-5,6] x [a</a,T\: \x\ < 



h 



(3.116) 



where we will choose h later. Note that h < const a for some (small) constant would 
assure S C T>(k). We will not impose such a restrictive condition on h but replace S by 
S H T>(k) in case S is not a subset of V(k). The following proposition gives our estimate 
on the first exit time from S. 

Proposition 3.14. Let t 2 c\yfo and (x2,t 2 ) G S. Then there exists a constant L > 
such that for any fx > ; we have 



F t2 > X2 {rs > t} < 
under the condition 



exp 



a(i, £2) 
1 + M e 



1 - 



1 



fJ>log(h/a) 



(3.117) 



(3.118) 



Proof: The proof follows along the lines of the one of [BG, Proposition 4.7], the main 
difference being the quadratic behaviour a_i 2 ?J a(t) ^ a + t 2 of a in our case as opposed 
to the linear one in [BG|. 

We start by defining a partition t 2 = uo < ■ • • < wi? = i of [£2, given by 



e h 2 

a(iik,Uk-i) = (1 + /x)-log-2. 



for 1 < k < K 



2a(t,t 2 ) 



;i + /i)£log(/i 2 /cj 2 ) 



(3.119) 



On each interval [uk,Uk+i], we consider a Gaussian approximation (x t )te[u k ,u k +i] °f x *> 
defined by 



(3.120) 



where = Wt — W Uk - If \x s \\fa{s) ^ /i for all s G [u^itfc+i], then by ([TJ) and (| 
there is a constant M > such that 



I ry, I <" 



l.9o 



(x u ,u)2; u |e a(s ' u)/e du 



< Af 



/r 3 



a(tt fc ) 3 / 2 a(wjfe) 



a a(« fc+1 ,u fe )/e 



(3.121) 
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for all s S [uk> u k+l]i provided the condition 

„2 



a(uk) 
M V a(u k+ i) 



-a(u k+1 ,u k )/e + 4 



holds for all k. Now, 

a{u k ) 



ca+ a(u k+ i,u k ) / a + a(ufc+i, it*,) 

^1+ Q .Q I -I- I O ~ 



'2 



t 3 



(3.122) 



(3.123) 



where c is a constant satisfying ^ o!{t) ^ ct for all t G [0, T]. This shows that there 
exists a constant L > such that the condition ( |3.122| ) is satisfied whenever 



(- 



1 + 



f 3 

'2 



r 2 



(7 



2 ' 



(3.124) 



This condition is equivalent to (3.118). 

Assume \x Uk \yJ a{u k ) ^ h for the moment. Then, 

P^*«J sup |x a |v^W<4 ^P Ufc '^{|xW +i | V ^K^)^2/ l } 



4/;. 



2vr^ ) +1 a('Ufc + i) 



(3.125) 



where i4t +1 denotes the variance of x^ +1 . By partial integration, we find 



„(fc) 



(7 



e 2a(« fe+1 , S )/ £ ds ^ 

u k a(u k+ i) L 



3 2a(u fc+1 ,« fe )/e 



(3.126) 



Now, the Markov property yields 

^ X2 {t s >t}= F t2 ^\ sup \x a \y/K(s) ^ h) < TT 



fc=0 



and the bound ( 3.117 ) follows by a straightforward calculation. 



(3.127) 
□ 



The preceding proposition shows that a path starting in S is likely to leave S after 
a short time. We want to show that such a path (or any path starting in T>{n) \ S) 
is also likely to leave T>(k). For this purpose, we will again compare xt to a Gaussian 
approximation, given by 



dx" 



-a (t)x^dt + -^=dW t , 

e V 6 



(3.128) 



where ao(t) = na(t), so that f(x,t)/x ^ ao(^) m T>(k). Assume that x t2 > 0. Then 
x s ^ x® holds as long as x s neither leaves T>(k) nor crosses the t-axis, cf. [ BG| , Lemma 4.8]. 
Therefore we can proceed as follows. Once a path is in T>{k) \S, there are two possibilities. 
Either, x® does not return to zero, or it does. If X® does not return to zero, then it is 
likely to leave T>(k) via the upper boundary and so is xt- So we are left with the case of 



36 



x® returning to zero. This event has a small but not negligible probability. Note that if 
returns to zero, then xt is still non-negative. If xt has nevertheless left T>(n), we are done. 
If not, xt is either in S or in T>(k) \ S. Since we may assume that, after a short time, xt 
is in T>(k) \ S again, we can repeat the above argument. 

Making the above-said precise, we obtain an integral equation for an upper bound on 
the probability that x u does not leave T>(k) up to time t, which will be solved by iterations. 
We will cite the integral equation from |BG| , as the general arguments leading to it do not 



require adaptation. Let us first introduce the necessary notations. We choose h = Ka for 
some (possibly large) constant K > 0. For k £ (0, 1), we choose n > in such a way that 



1 ii 



21 + n 



1 - O 



1 



fi logK 



1 - o 



1 



//log-fT 



(3.129) 



for all large enough K. Note that choosing k too close to 1 requires large [i and is thus 
not desirable. Since we want to apply Proposition |3. 14 on the first exit time from S with 
t 2 ^ C2\/a for a suitably chosen c 2 , Condition ( 3.118 ) must be satisfied. Therefore, we 
choose c 2 = C2(K) large enough for 



to hold. Now, set 



g(t,s) 



-Ka(t,s)/e 



-2Ka(t,s)/e 



and 



C 



maxi 



rx(t)V^) A and 
I \/ira > 



,(o) 



7TKI — 



-Kh 2 /a 2 



(3.130) 



(3.131) 



(3.132) 



For t 2 s$ s ^ t < T, let Ql v '{s) = 1, and define Qj n) (s) for n ^ 1 by 



Q[ n \s) = Cg(t, s) + ce - K «(*' s )/ £ + c f Q^" 1 )^)^! e -™M/e ^ + g ( Ui s )j dn _ (3 _ 133) 



Then, [BG, (4.95) and (4.107)] show that for any n and s ^ t 2 ^ c 2 yfa, 



sup F s > x {t v{k) >t}^2 

x: (x,s)s5 



a 



-KOt(t,s) / E 



and 



sup J 

x: (x,s)£V(k)\S 



Q 



~'{Tv( K )>t}^Ql n) (s) 



(n) 

Next we estimate Q t by showing that 

,(«) 



Q\ n} (s) ^ Cg(t, s) + a n e ~ Ka ^/ 2£ +b n V 



n 



holds with ai = c, b\ = 3c/ k in the case n = 1, and with 



/3c 



4C\ 



n-2 



6c\i /6c\ n 



-1 



j=0 



4C 



c / 4C 

^^ 2C ( 1 + -, 



(3.134) 
(3.135) 

(3.136) 

(3.137) 
(3.138) 



37 



for n > 1, provided 6c/k ^ 1/2. Note that the latter imposes a condition on K = h/a. 
To obtain the bound ( 3,136| ), we proceed as in |BG|, the only difference lying in the term 
a n e~ Ka ( t,s ^ 2e , where we sacrifice a factor of 2 in the exponent in order to gain a smaller 
coefficient a n . Our choice of a n yields a less restrictive condition on h/a, namely we only 
need 



12^ shy+» c _ Kh 2 /a 2 ^ l 



(3.139) 



which is satisfied whenever K = h/a is large enough. 



Now, ( 3.137 ) and ( 3.138| ) imply that for K and c<z(K) large enough, 

for all t ^ t2 ^ c^^fa 

(3.140) 



sup ^{r v{K) > t} ^ C (-=) -= 



rea(t,t 2 )/e 



with some constant Cq. This completes our outline of the proof of Theorem |2.4|. 



4 Asymmetric case 

We consider in this section the nonlinear SDE 



dx t = -f(x t ,t)dt + ^ F dWt, (4.1) 

e ye 



where / satisfies the hypotheses given at the beginning of Subsection 2.2. By rescaling x, 
we can arrange for d xx f(x c , 0) = —2, so that Taylor's formula allows us to write 

d x f(x c + x,t) = d x f(x c ,t) + x[-2 + n(x,t)] (4.2) 

where n € C 1 and n(0, 0) = 0. Since d x f(x c ,t) = 0(t 2 ) by assumption, d x f(x,t) vanishes 
on a curve x (t) =x c + 0(t 2 ). We further obtain that 

f(x(t) + z,t)= f(x(t),t) + z 2 [-1 + r (z, t)] 

f(x(t),t) = f(x c , t) + 0(t 4 ) = a + ai t 2 + 0(t 3 ), 

where tq £ C 1 and ro(0, 0) = 0. Thus f(x, t) vanishes on two curves x*,(t) and Xg(t), which 
behave near t = like x c ± V^o + ait 2 [l + 0(vV+ai?)], as indicated in fl2.34| ). The 
behaviour of the linearization follows from (|4l 

4.1 Deterministic case 



The proof of Theorem 2.5 follows closely the proof of Theorem |2.1| , with some minor 
differences we comment on here. The dynamics of yt = %t ~ x +(t) is still governed by an 
equation of the form 

e^ = aX(t)y + bX(y,t)-e^±, (4.4) 
but now Taylor's formula yields the relations 

<(t)x-(V^V|i|) (4.5) 

b%(y, t) = -y 2 [1 + O(v^) + 0(t) + 0(y)] , (4.6) 

while relation ( 3.10| ) holds for the derivative of x* + , with t* replaced by ti_. Lemma |3 
becomes 
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Lemma 4.1. Let a(t) be a continuous function satisfying a(t) x — (/3 V \t\) for \t\ ^ T, 
where (3 = (3(e) ^ 0. Let \o ~ 1> anc ^ define a(t,s) = a(u)du. Then 



<5(t-T)/ e+ i /"* e <5(«,s)/e ds 



i 



1 



/or \t\ </3V Vi 

for (3\J ^/e^ \t\ ^ T. 



(4.7) 



Proposition 3^3 carries over with some obvious adjustments, and shows the existence 
of a constant cq such that 



Vt 



for — T ^ t < t = -co(a/oo V a/e). 



(4.8) 



In particular, yt x (e/y/ao) A \/e. An adaptation of Proposition 3/4 yields the existence 
of a constant 70 > such that, for ao ^ 7q£, 



lfe = Ci(t)(t^.-t) + C72(t) with Ci(i)x-, C 2 (t)xi- 



(4.9) 



for all |t| ^ I to |- This shows in particular that yt vanishes at a time t satisfying t — t* 
e/ yfaQ. Proposition [T5] is replaced by 



Proposition 4.2. Assume that a$ < jo£. Then, for any fixed t\ x yfe, 

x t - x c x yfe fort ^t^ t 1} 
and xt crosses x*\_(t) at a time t satisfying t x yfe. 

Proof: Let x t = xt — x c . We first observe that by Taylor's formula, 

f{x c + x,t) = f(x c , t) + xd x f(x c , t) + x 2 [-1 + O(x) + 0(t)] . 
This shows that 

£ ^ = a + ai t 2 -x 2 + 0(t 3 ) + 0{t 2 x) + 0{ti 2 ) + 0(x 3 ). 
Thus, with the rescaling 

1/4 r , -1/4 r 

x = a{ y/ez, t = a 1 yes, 



we obtain that zt obeys a perturbation of order yfe of the Riccati equation 

with oq 



dz 22 
— = a + s - z , 
ds 



1 Qq 7o 



(4.10) 



(4.11) 



(4.12) 



(4.13) 



(4.14) 



One easily shows that the solution satisfies z s x 1 for s of order 1, and this property 
carries over to the perturbed equation with the help of GronwalFs inequality. Finally, since 
Xt ~ an£ i x +(t) ~ x c x \/flo V |*|, these curves necessarily cross at a time i X \/e. □ 
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The assertion on the existence of a particular solution x tracking the unstable equi- 
librium branch Xg(t) follows from the observation that z s = x_ s satisfies the equation 

dz 

e^ = -f(z s ,-s). (4.15) 

This system admits Zq(s) = Xq(—s) as a stable equilibrium branch. Thus the same ar- 
guments as above can be used to show the existence of a solution z s tracking Zq(s), with 
similar properties. Proposition |3.7| admits the following counterpart: 

Proposition 4.3. For all t G [— T,T] and all ao = O e (l), 

a(t) := d x f(xf et , t) x -(|t| V^Vv^) (4.16) 
a(t):=d x f(xf c \t) x \t\ V^V y/i. (4.17) 



Proof: The proof is a direct consequence of ( |4.2| ) and the properties of xf et , and thus 
much simpler than the proof of Proposition |3.7j. □ 



Finally, with Lemma |4.l| , we immediately obtain 



2\a( 



3 2a(t,a)/e ^ g 



1 



(4.18) 



4.2 The transition regime 

We consider now the regime of a sufficiently large to allow for transitions from the potential 
well at x\ to the potential well at x* , by passing over the saddle at Xq. Here x^ et and 
Xi et denote solutions of the deterministic equation 

e^ = f(x,t) (4.19) 

tracking, respectively, the stable equilibrium branch and the unstable equilibrium 

branch x^(t), while x% denotes a general solution of the SDE ([4.1|). Our aim is to establish 
an upper bound for the probability not to reach a level 5q between x^(t) and x* (t), situated 
at a distance of order 1 from both equilibria. [ |BG| , Theorem 2.3] shows that if xt reaches 
5q, and 5o is close enough to x*_{t) (but it may still be at a distance of order 1), then it is 
likely to reach a neighbourhood of x* (t) as well. 

Let 5q < $i < x c < 82 be the constants satisfying ( 2.46Q , that is, 



f(x,t) X —1 for Jn ^ x ^ 61 and \t\ ^ T 

d X xf{x,t) ^ for 5\ ^ x ^ 62 and |i| ^ T. 



The basic ingredient of our estimate is the following analogue of Lemma 3.11 



Lemma 4.4. Fix some initial time to G [— T, T\. We consider the following processes on 
[to,T]: 

• the solution xf et of the deterministic differential equation fl4.19| ) with an initial condi- 
tion x^f G [61,62], such that xf ct ^ Si for all t G [t ,T]; 

• the solution xt of the SDE ( f4.1| ) with an initial condition xt G [xf ct , 62]; 

• the difference yt = %t — xf et , which satisfies y to = xt — xf^ 0; 
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the solution y® of the linear SDE 



dy° = ±a(t)y°dt + ^dW t , 
e ye 



where a(t) =d x f(xf et ,t) 



(4.21) 



„detl 



with initial condition yf Q E [yt , 82 — x± j . 

If <5i ^ Us +Xg Ct ^ 82 for all s E [to,t]> 2/s ^ 2/s / or ^ose s. Similarly, if 8\ ^ x s ^ ^2 
/or aZZ s E [iojt]> ^en y|? ^ y s /or i/iose s. The result remains true when t is replaced by 
a stopping time. 

We will proceed as follows. Let £(t) be the function denned in ( [4.18 ), and let h be 
such that Xg et + h\J £(s) < 62 for all s E [to?*]- Given x E (81,82) and times to < ti < t 
in [— T, T], we consider the solution xt of the SDE ( f4.1[ ) with initial condition xt = x$. 
We introduce the stopping time 



r = inf{s E [to,ti]: x s < 5i} E [*o, *i] U {°o}- 



(4.22) 



We can thus write 



¥ t0 ' X0 {x s > 5q Vs E [t ,*]} < P* ' 10 !^ > 5! Vs E [t ,ti]} 

+ E^{ W ^{x. >5 Vs E [r,t]}}. (4.23) 

The first term can be further estimated by 



>{x s ><5i VsE [t Q ,ti]} <P'°'*°{ sup 



det » 

- > h\ 



+ F t °' x °{8 1 < x s < 4 et + /i v / CM Vs E Mi]}- (4.24) 

The two summands in (|4.24j) can be estimated in a similar way as in the symmetric case. 
The first one is dealt with in the following result. 

Proposition 4.5. Assume 6% ^ xq ^ xf et + ^/i\/C(*o)- T/ien 



p*0,£0 I 



„dct 



sup 



> h 



}< ( 



5 /|a(ti,t )| 



+ 1 ) e-^> 2 



(4.25) 



where k is a positive constant and a(ti,to) = f t * a(s) ds. 



The proof is almost the same as the proof of Proposition 3.12. Instead of ( 3.84| ), we 
may define r° and r + as the first times when ys either reaches Si — x^ et or the upper 
boundary h-\/((s). Then Lemma L4 implies y s ^ yi for s E u^+i Ar°Ar + ]. However, 
when estimating the probability that r° < u^+i as in ( 3.88p , it is sufficient to use the fact 
that T° is larger than the first time y^ reaches 0. Finally, ( |3.91| ) still holds with the 
present definitions of £ and a, because of (4.18]) and Proposition [47 
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Let us now examine the second term in ( 4.24 ). 



Proposition 4.6. There exist constants c\ > and k > such that, if c^ 2 o~ a^ 4 Ve 3 / 4 
and > 2<t, then 



• xo {dt <x s ^ x d s ct + h^/(^) Vs g [-do^M]} 

<-ex ( R 1 a(ti,-ci(7 2 / 3 ) - 

^2 6XP \ K log(/i/V) V |log<r| e 



(4.26) 



/jo/ds with a(t,s) = J s a(u)du, for —c\o~ 2 ^ ^ t\ ^ cio -2 / 3 and aZZ initial conditions xq 
satisfying 5i < x ^ ^^2/3 + /i\/C(-eio- 2 / 3 ). 

Proof: 

1. Let a? det be the deterministic solution tracking the saddle at Xg(t) and set xt = xf et +z t . 



Then 



1 



dz t = - [a(t)z t +b(z t ,t)] dt + dW t , 



where ( p7|) , flOp and (pL2C|) imply 



a(t) x |t| V ^/ao V 



at) 



and 



6(z,t)s^0 for xf ct + ,z G [tf l9 tf 2 ] 



(4.27) 



(4.28) 



(4.29) 



Let £ = g(h/a) ^ 1 and define a partition — ci<7 2 / 3 = no < • • • < n# = ii of [— ci<r 2//3 , i] 

by 



a(uk, Uk-i) = ge for 1 ^ k < K 



a(ti, -ci<r 2//3 ) 



Q£ 



(4.30) 



Let 
Qk '- 



sup 



< x det + z s ^ * det + /iy 7 ^) Vs G [n fc ,n fc+1 ]}. 



<5l <2 dot +Zk ^det -ffty^TJ 

Then we have, as in ( 3.99| ), 

p-^'^fr <x s ^ x dct + h^/Os) Vs G [-c^M] ^l\Q k 



K-l 



(4.31) 



(4.32) 



k=0 



The result will thus be proved if we manage to choose g in such a way that Qk is 
bounded away from 1 for k = 0, . . . , K — 2. 
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2. We will estimate the Q k in a similar way as in Proposition 3.13 , but we shall distinguish 
three cases instead of two. These cases correspond to x s reaching the levels x det , x det 
and 5\. We introduce a subdivision u k < u k \ < u kj2 < u k+ \ defined by 



1 



a(uk,2,Uk) 



and stopping times 

r fci i = inf{s G [u k ,u k ,i] ■ z s ^ x det 
Tk,2 = inf{s G [u k ,u kj2 \- z s ^ 0}. 

Then we can write, similarly as in ( 3.1 10|) , 



>dct1 



t> u k, z k 



{5 1 < x dct + z s ^ x det + h^/Os) Vs G [u k ,u k+1 ]} 



(4.33) 



(4.34) 



(4.35) 



< p« fc A-{ x det < jdet + ^ < x dct + ^y^y Vg G [ Uk} u k>1 ]} 

+ E u ^{l {Tki<ilk i} F T ^ z ^{5 1 < x det + z s ^ x** + h^/Os) Vs G [r kA ,u k+1 ]}} 



The first term can be bounded by comparing with the solution of the SDE (4.1) 
linearized around xf et , with the help of Lemma 4.4. As in (3.112), we obtain the 
upper bound 



V 1 



h 



-2q/3L 



sup 



V|a(n)|CK)-e-^ 3i , 



a 



(4.36) 



where L > is a constant such that a(u) ^ L|a(u)| for all u. Now if r k> i < u ki \, we 
also have 



^'^{Sx < x det + z s ^ x d s et + hy^Of) Vs G [r kjl ,u k+1 ]} 
<: {o < z s ^ x det - x dct + h^Osj Vs G [r M , 



(4.37) 



E 



{Tfc,2<«fc,2} 



I)T"fc,2,2T,. 



< x dot +z s < x det + /i V / C(^) Vs G [r fc , 2 ,n fc+1 ]}}. 



Comparing with the solution of the SDE ( |4.1| ) linearized around x det , the first term 
can be bounded, as in (3.104), by 



1 



7T V 1 



-2£>/3 



sup 



V a ( n ) / det ~det\ 

sup (x u -x u ). 



a 



(4.38) 



This estimate shows that a path starting on x dct at time T k i has an appreciable 
probability to reach the saddle before time u k 2- Note, however, that we cannot obtain 
directly a similar estimate for the probability to reach 8\ as well, which is why we 
restart the process in x det . 
3. In order to estimate the second summand in ( [4.37| ) , let be the process starting in 
at time r k 2 which satisfies the linear equation 



dz^ = -a(s)zW ds + 4= dVFj fc \ 
e ye 



(4.39) 
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with = W S — W Tk 2 . The variance VuJ +1 of z^J +1 satisfies, as in ( |3.101[ ), 



(fc) 



.(*) 



e u u k+1 & 



1 - e" 2 ^ 3 c a 2 
mi 



« fe ^«^« fe+ l a(li) 



(4.40) 



Thus we obtain, using Lemma 4.4, 



1 



fftdet _X,-\["„W "1-1/2 



11/1 

^2 + 7?vr^ 



-V3 



sup 



-{x 



"fe+1 



<5i) e 



-e/3 



(4.41) 



Here the introduction of the stopping time 2 turns out to play a crucial role. The 
above probability is indeed close to 1/2 when £ is larger than a constant times | log <r | , 
which shows that once a path has reached the saddle, it also has about fifty percent 
chance to reach the level Si in a time of order e|log a\/a(u). 
4. From (pUffi ), (JO|) and (|4.41|) and the fact that g > 1, we obtain the existence of a 
constant C > such that 



Qk < 7T + 



2 CT u k f^uf^u 



fc+i 



(4.42) 



+ sup (x* 



Since |u| cio -2 / 3 , the properties of x det , x det , £ and a imply the existence of another 
constant C\ such that 

e -e/3- 



\ + Ci 



• • ,2/3 J" ^ 

Now choosing cf 2 = l/(18Ci) and g = 3£log(18Ci/i/cr) V 31og(18Ci^/cI/a 2 / 3 ) V 1, 
we get Qk ^ 2/3 for k = 0, . . . , K — 2 and thus 



fc=0 



2 (3/2) 



A" 



^ - exp< 



log- 



3^ q(ti,-ci(7 2 / 3 ) i 

Q£ J ' 



and the result follows from our choice of g. 



(4.44) 



□ 



It remains to estimate the second term in ( 4.23 ), describing the probability not to 
reach So when starting in 5\. This is done by using the fact that, by assumption, the drift 
term is bounded away from zero on the interval [So, S\]. We will need to assume that it can 
be bounded away from zero on a slightly larger interval, which is possible by continuity of 
/• 

Proposition 4.7. Let < p ^ Si — So be a constant such that f(x,t) ^ —fo < for 
So ^ x ^ Si + p and \t\ ^ T. Then 

P* ' 51 {x s > So Vs G [t , to + ce] } < e~ k/a2 (4.45) 
holds for all to £ [—T,T — ce], where k = fop 2 /4(Si — So) and c = 2(Si — So)/ fo- 
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Proof: Let x^ be defined by 



x° = 5 1 -—(t-t ) + -^=W t - to , t^to. (4.46) 
e \ e 



By Gronwall's inequality, it is easy to see, as in Lemma i.4 , that if 5q ^ x® ^ 5\ + p for 
all t E [to, to + ce]j then x t ^ x° for those t. We thus have 



F t0 ^{x s >6 Vs e [t , to + ce]} < P* oA { sup x° s + ^(s - t ) > & + p) 

+ P* oA {5 <x°<5i + p- j(s-to) Vs 6 [t ,t + ce]}. (4.47) 
Note, however, that for s = to + ce, 

& + P ~ j(s - to) = St + p - 2(5 l - So) < <y 0j (4.48) 



so that the second term in ( 4.47| ) is equal to zero. The first term equals 



dO.O 



by Doob's submartingale inequality. □ 
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